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Abstract 


This  report  Is  concerned  with  the  measurement  of  frequency  (Doppler 
shift)  by  a collection  of  sensors  which  may  be  grouped  into  two  or  more 
clusters  (subarrays).  The  signal  may  be  a sinusoid,  or  it  may  be  a narrow- 
band  Gaussian  process  whose  center  frequency  is  then  the  quantity  of 
interest.  For  much,  chough  not  all,  of  the  analysis,  it  is  assumed  that 
Che  noise  is  Incoherent  from  sensor  to  sensor.  The  signal  wavefront  is 
assumed  to  be  perfectly  coherent  over  all  sensors,  but  the  angle  subtended 
at  the  source  by  the  various  sensors  may  be  sufficiently  large  so  chat 
differencial  Doppler  shifts  cannot  be  Ignored.  The  following  are  some  of 
Che  more  general  conclusions. 

I.  No  differential  Doppler  shift  between  sensors 


1. 


Sinusoidal  signal 

-3 

a.  The  mean  square  frequency  estimation  error  varies  as  T , where  T 
is  the  observation  time.  It  varies  as  the  Inverse  second  power  of 
Che  posc-beamforming  signal  to  noise  ratio. 

b.  The  frequency  estimate  is  not  degraded  by  lack  of  prior  knowledge 
of  bearing  and  range  (i.e.  sensor  Co  sensor  delay  of  the  signal). 
Frequency  estimation  and  bearing-range  estimation  are  therefore 
uncoupled. 

c.  The  aaan  square  astlmatlon  error  of  sensor  to  sensor  signal  delay 
varies  as  T~^.  Hence  the  frequency  estimate  becomes  extremely 
•ecurate  long  before  the  delay  estimate. 


2.  Narrowband  Gaussian  signal,  TW  >>  1. 


The  frequency  of  Interest  Is  here  the  center  frequency  of  the  narrowband 
Gaussian  signal.  This  particular  computation  has  been  carried  out  for 
completely  arbitrary  array  geometry  and  noise  field  properties. 

a.  The  mean  square  estimation  error  of  the  center  frequency  varies 
as  T Accuracy  therefore  Improves  at  the  same  time  rate  as  for 
delay  estimation  but  much  slower  than  for  frequency  estimation 
using  a pure  sinusoid. 

b.  The  mean  square  estimation  error  of  Uq  depends  on  two  quantities 
aside  from  T:  the  post-beamformlng  signal  to  noise  ratio  of  the 
beamformer  optimally  matched  to  the  noise  field  and  the  slope  (on 
a db  scale)  of  the  signal  spectrum. 

c.  The  coupling  between  center  frequency  estimation  and  bearing 
estimation  depends  on  the  derivative  of  the  array  gain  with  respect 
to  bearing  In  the  neighborhood  of  the  true  bearing.  If  this  deriva- 
tive Is  zero  there  Is  no  coupling.  Thus  bearing  and  frequency 

estimation  are  completely  uncoupled  for  a spatially  Incoherent  noise  ^ 

field.  They  are  only  weakly  coupled  In  many.  If  not  most,  practically 
Interesting  situations.  For  strong  coupling  the  noise  field  would 
have  to  exhibit  a rapidly  varying  spatial  structure  near  the  signal 
bearing  (e.g.  an  Interference  close  to  the  signal  In  bearing). 

3.  Narrowband  Gaussian  signal,  TW  « 1. 

When  the  tlma-bandwldth  product  is  somII  compared  to  unity,  a Gaussian 
signal  looks  essentially  like  a sinusoidal  signal  of  unknown  amplitude 
and  phase.  The  amplitude  estimate  Is  not  coupled  to  the  frequency 
estimate  so  that  the  results  of  section  (1)  are  applicable.  One  must 


note,  however,  that  the  apparently  very  accurate  frequency  estimate 
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(T  dependent)  is  a short  term  frequency  estimate.  To  define  a center 
frequency  Uq  of  the  narrowband  signal  one  Inherently  requires  an 
observation  time  large  compared  with  the  Inverse  signal  bandwidth. 


II.  Differential  Doppler  estimation. 


The  differential  Doppler  estimation  problem  was  studied  for  two  or  three 
sensor  groups  (subarrays)  which  were  treated  as  equivalent  single  sensors. 

1.  Sinusoidal  signal 

-3 

The  mean  square  error  for  differential  Doppler  estimation  varies  as  T , 
just  as  that  of  the  frequency  estimate.  For  two  subarrays  It  differs 
from  the  latter  by  3db  because  there  are  two  opportunities  to  measure 
frequency  but  only  one  to  measure  frequency  difference.  Subtraction  of 
two  separate  frequency  estimates  to  obtain  the  differential  Doppler 
shift  leads  to  the  same  error  as  direct  estimation  of  the  differential 
Doppler  shift. 

2.  Narrowband  Gaussian  signal. 

Since  narrowband  signals  with  TW  <<  1 have  the  properties  of  sinusoids, 
only  the  case  TW  >>  1 was  considered.  An  analytically  Important  assump- 
tion places  an  upper  bound  on  values  of  T for  which  the  results  are 
valid:  If  the  signal  bandwidth  Is  o,  the  differential  source  velocity 
Av  and  tha  velocity  cf  sound  c,  T must  be  shorter  than  the  period  of  a 
sinusoid  of  frequency  (Av/c)o. 

a.  The  mean  square  estimation  error  of  differential  Doppler  shift 
-3 

varies  as  T , Just  as  the  frequency  estimate  of  a sinusoid.  The 
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p 
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estimation  process  for  differential  Doppler  shift  Is  coherent 
whereas  that  for  center  frequency  Is  Incoherent. 

b.  Vfhen  the  signal  to  noise  ratio  In  the  signal  band  Is  high,  the 
differential  Doppler  estimate  using  narrowband  signals  has  the 
same  mean  square  error  as  an  estimate  using  sinusoidal  signals. 

c.  The  only  Important  difference  between  differential  Doppler  estimates 

using  narrowband  signals  at  high  and  low  signal  to  noise  ratios  Is 

the  slgnal-to-nolse  (S/N)  dependence  of  the  mean  square  error.  At 

high  signal  to  noise  ratios  It  varies  as  (S/N)  at  low  signal  to 

_2 

noise  ratios  as  (S/N) 

d.  Because  of  the  different  T dependence,  efforts  to  obtain  differential 
Doppler  shift  from  separate  measurements  of  center  frequency  at  each 
subarray  lead  to  a much  poorer  estimate.  The  loss  In  performance 

Is  proportional  to  the  square  of  the  TW  product. 

e.  The  estimate  of  differential  Doppler  shift  Is  completely  uncoupled 
from  the  estimate  of  bearing,  range  and  parameters  describing 
spectral  shape.  Differential  Doppler  measurements  can  therefore  be 
carried  out  without  concern  about  most  other  parameters  which  are 
likely  to  be  unknown  In  practice. 

f.  When  there  are  three  subarrays,  there  are  two  differential  Doppler 
shifts.  Their  estimates  are  coupled,  but  neither  Is  coupled  to  the 
estimates  of  bearing,  range  or  spectral  parameters.  When  the  signal 
to  noise  ratio  In  the  signal  band  Is  high,  the  estimation  error  for 
each  differential  Doppler  shift  is  the  same  as  that  from  the  appropriate 
pair  of  subarrays.  When  the  in-band  signal  to  noise  ratio  is  low,  a 
small  gain  can  be  ude  by  processing  the  three  subarray  outputs  as  a 
unit  rather  than  obtaining  separate  differential  Doppler  shifts  from 
two  pairs  of  subarrays. 


i 
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I.  Introduction 


This  report  sumnarlzes  work  carried  out  under  contract  N66001-76-C-0182 
between  Yale  University  and  the  Naval  Ocean  Systems  Center,  San  Diego.  It 
describes  the  continuation  of  studies  Initiated  under  contract  N66001-75-C-0210 
and  reported  In  June  1976.  [1] 

The  overall  objective  of  both  contracts  was  to  obtain  Information  about 
an  acoustic  source  by  means  of  a large  array  of  sensors  which  may  be  grouped 
Into  several  subarrays.  Reference  [1]  was  concerned  with  source  localization 
In  bearing  and  range.  Here  we  turn  to  the  problem  of  Doppler  estimation  and 
Its  coupling  to  bearing  and  range  estimation. 

Much  of  the  basic  theory  developed  In  [1]  remains  applicable  to  the 
present  study  and  will  therefore  only  be  referenced  briefly  when  needed. 

In  particular,  we  shall  retain  most  of  the  basic  assumptions: 

(a)  The  source  radiates  either  a sinusoid  or  a narrowband  Gaussian 
signal. 

(b)  The  signal  wavefront  Is  perfectly  coherent  over  the  entire 
receiving  array. 

(c)  The  noise  Is  Gaussian  with  a bandwidth  broad  compared  with  that 
of  the  signal.  Much,  though  not  all,  of  the  analysis  will  also 
assume  that  the  noise  Is  spatially  Incoherent. 

(d)  The  frequency  or  center  frequency  of  the  signal  observed  at  any 
given  sensor  does  not  change  significantly  over  the  observation 
period. 

However,  In  contrast  to  the  earlier  work,  the  signal  frequency  or 
center  frequency  need  not  be  the  same  at  each  sensor.  Thus  we  allow  Doppler 
shifts  and  even  differential  Doppler  shifts  from  sensor  to  sensor.  Assump- 
tion (d)  does  not  preclude  source  motion;  it  merely  says  that  the 
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source-receiver  geometry  does  not  change  significantly  during  the  obser- 
vation period. 

A major  conclusion  of  the  previous  report  was  that,  for  most  purposes, 
each  subarray  could  be  treated  as  a single  sensor  whose  relatively  weak 
directional  properties  were  of  only  secondary  importance.  Thus  one  could, 
without  significant  loss,  beamform  at  each  subarray  and  treat  the  resulting 
beams  as  the  outputs  of  equivalent  single  sensors.  The  subarrays  serve  only 
to  enhance  the  signal  to  noise  ratio  of  these  equivalent  sensors.  If  the 
source  is  sufficiently  distant  to  warrant  the  use  of  more  than  one  subarray, 
there  will  be  very  little  differential  Doppler  shift  over  the  elements  of 
any  one  subarray.  Beamforming  at  the  subarray  level  should  therefore  not 
destroy  any  significant  amount  of  Doppler  information.  We  shall  find  this 
to  be  true  in  the  first,  simple  formulations  of  the  Doppler  estimation 
problem.  Thereafter  we  shall  then  treat  each  subarray  as  a single  equivalent 
sensor  with  appropriately  enhanced  signal  to  noise  ratlov 

As  in  Kef.  [1]  we  shall  use  the  Cramer-Rao  inequality  to  set  lower 
bounds  on  attainable  estimation  errors.  In  the  limit  of  large  observation 
times,  this  lower  bound  is  approached  by  the  error  of  a maximum  likelihood 
estimator  [2].  At  least  in  this  important  limit,  therefore,  the  calculated 
bound  can  be  reached  by  an  actual  instrumentation. 


In  the  most  general  setting 
of  the  problem  we  have  a geometry 
such  as  shown  in  Figure  1. 

A source  located  at  range  r 
and  bearing  0 (measured  in  some 
convenient  polar  coordinate  system) 
radiates  a signal  which  is  received 
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(together  with  the  Inevitable  noise)  at  M sensors  in  known  but  perfectly 
arbitrary  locations.  The  output  of  the  1^^  sensor  is 

XfCt)  - 8^(t)  + n^(t)  (] 

s^(t}  is  the  signal  component  and  n^(t)  the  noise  component.  We  construct 
a data  vector  x as  follows: 

(1)  Represent  x^(t)  by  a set  of  n numbers  [e.g.  the  time  samples 
Xi(ti),  x^(t2) . . . . iX^(t^)  or  the  Fourier  coefficients  c^^^, 

‘'2i’-**’‘'ni^- 

(2)  Concatenate  these  n number  sets  for  all  index  values  i [e.g. 

Xl(ti)  ’ • • • ’ ^2^*^!^ . • • • .X2  (tjj) , . . -XjjCtj^) , . . .Xjj(t^)  ] . 

The  resultant  data  vector  x therefore  has  dimension  nM. 

The  Cramer-Rao  technique  requires  computation  of  the  likelihood  function 

A - p(x/a^,a2, . . .0^)  (2) 

02,... >a^  are  the  parameters  to  be  estimated  (e.g.  Doppler  shift,  bearing, 
range,  etc.).  The  right  side  of  Eq.  (2)  is  the  conditional  probability 
density  for  the  data  vector  x when  the  unknown  parameters  have  values 

Ul»®2 * • • *'*1'  * 

Next  one  must  form  the  Fisher  information  matrix  J whose  elements  ^ 
are  given  by 

E{  } is  the  expectation  of  the  bracketed  quantity.  The  Cramer-Rao  In- 

2 ^ 

equality  asserts  that  the  (linbiased)  mean  square  estimation  error  D (o^) 
of  the  parameter  is  bounded  by  the  kk  element  of  the  Inverse  matrix  J 

i (4) 

Off-dlagnol  elemanta  of  J ^ measure  the  coupling  between  estimation  errors. 
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II.  Frequency  Estimation  - No  Differential  Doppler  Shift 

We  begin  with  the  problem  of  frequency  estimation  using  a single  sub- 
array. The  elements  of  such  an  array  are  presumably  spaced  rather  tightly 
and  differential  Doppler  shifts  from  sensor  to  sensor  are  negligible  for  all 
but  very  nearby  sources.  The  frequency  or  center  frequency  of  the  signal 
component  received  by  each  sensor  is  therefore  the  same.  Our  objective  is 
to  estimate  that  frequency  or  center  frequency,  possible  in  the  presence  of 
other  initially  unknown  parameters  such  as  source  bearing  or  range. ^ 

1.  Sinusoidal  signal 

We  begin  with  the  simplest  case,  a source  radiating  a pure  sinusoid. 

The  Doppler  shifted  signal  component  at  the  receiving  array  is  a sinusoid 
of  frequency  Thus  the  output  of  the  i sensor  [Fig.  1]  is 

x^(t)  - sin[a)^(t-T^)  - ((i]  + n^(t) 

is  the  signal  travel  time  from  the  source  to  the  i^^  sensor  (bearing  and 
range  dependent)  and  n^(t)  is  the  noise  component  at  the  i^^  sensor.  For 
the  present  computation  we  shall  assume  that  the  noise  is  Gaussian^  spatially 
incoherent  (statistically  independent  from  sensor  to  sensor)  and  spectrally 
white  over  a band  of  W Hertz  beginning  at  zero  frequency. 

Calculation  of  the  Fisher  information  matrix  proceeds  much  as  in 
Ref.  [1].  The  data  vector  consists  of  time  samples  of  the  sensor  outputs 
taken  at  the  Nyqulst  rate  (i.e.,  at  intervals  of  1/2W  seconds).  Because 
the  noise  is  both  spatially  and  spectrally  white  as  well  as  Gaussian, 

^Note  Chet  this  formulation  of  the  problem  does  not  permit  separate  esti- 
mation of  radiated  source  frequency  and  Doppler  shift.  Only  the  sum  of 
Che  two  frequencies  is  accessible  to  measurement  and  hence  to  estimation. 


(5) 
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ail  of  the  samples  are  statistically  independent  and  the  likelihood  function 
assumes  the  simple  form 


M n 1 2 

A - C exp  -{Z  Z -^iTrTrtx.Ct,)  - 8,(tJ]  } 


i-1  j-1  ^^i”  ^ J ^ 


(6) 


,th 


is  the  spectral  level  of  the  noise  received  at  the  i sensor  and  N^W 

is  the  average  noise  power  of  each  sample  at  that  sensor.  C is  a normalizing 

constant  which  insures  that  A is  a probability  density.  Note  that  J serves 

as  an  index  for  the  time  samples,  1 for  the  sensor  locations  (space  samples). 

As  a typical  element  of  the  Fisher  information  matrix  consider 
2 2 

-E{3  log  A/3u^  }.  Differentiating  the  natural  logarithm  of  Eq.  (6)  twice 
with  respect  to  one  obtains: 


M n 


38, (tj 


3 s,(t,) 


3w  2 N^W  3^,^  i j i j 2 


(7) 


,th 


,th 


x^(tj)  - s^(tj)  is  simply  the  J sample  of  the  noise  at  the  1 sensor. 
Since  the  noise  has  zero  mean,  the  averaging  operation  eliminates  the 
second  term  of  Eq.  (7)  and  one  is  left  with 

.2,  , M n , 3s, (t.)  ^ 


3uj 


i**l  J“1  3u)q 


(8) 


From  Eq.  (3) 


SfCt^)  - sinL^Ctj  - T^) 


(9) 


Hence 


2 1 » M 

a log  A 

-E{  ^ 2 ^ ^ “ N W'^J  *i 

1-1  J-1  1 ■’ 


3u 


n 

Z 


C08^Ia.^(tj-T^)  - 


Wlien  Che  noise  bandwidth  W is  sufficiently  large  so  that  there  are  many 
samples  in  a signal  period  one  can  approximate  the  J sum  in  Eq.  (10)  by 


(10) 


an  Integral. 
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M'  'A, 


_E^i — - 2W  -T.)^  cos^[u  (t  -T.)  - (>i]  At 

3w  ^ i-1  J ^ o J 1 


M A 
i-1  i 


2 t -^T 


2 2 

(t-T^)  cos  [u^(t-T^)  - <p]dt 


(11) 


In  the  first  version  of  this  equation  At  • 1/2W  is  the  time  increment  from 
sample  to  sample. (t^,  t^+T)  is  the  observation  interval. 

In  practice  T is  almost  certainly  very  large  compared  with  the  signal 
period.  In  that  case  double  frequency  terms  make  negligible  contributions 
to  the  integral  and  one  obtains 

“ I Z ^ {(T+t  - T - (t  -T  )^}  (12) 

.2  3..N^  o i oi 


Suppose  the  signal  travel  time  from  the  source  to  the  center  of  the 

array  is  One  can  absorb  that  quantity  into  with  the  definition 

t,  ■ t -T  . The  differential  delays  t.-t  are  almost  certainly  small  com- 
1 o o i o 

pared  with  I and  one  can  then  write  to  an  excellent  approximation 

_E^iil2i_A}  .i  j +T)^  - t (13) 

.2  3 j 1 N.  1 1 

dw  i“l  i 

o 

If  u)^  is  the  only  unknown  parameter,  Eq.  (13)  is  the  only  element  in 
the  Fisher  information  matrix  and  its  reciprocal  should  be  the  minimum 
mean  square  estimation  error.  This  is  a correct  but  highly  misleading 
statement.  Eq.  (13)  depends  strongly  on  the  largely  arbitrary  t^.  One 
can  apparently  reduce  the  estimation  error  drastically  by  merely  starting 
the  observations  laterl  There  is  a simple  explanation  for  this  paradox: 

We  are  making  the  completely  unrealistic  assumption  that  the  phase  angle 
^ is  known.  If  this  were  actually  true  we  could  use  the  phase  shift 
w^t^  accumulated  prior  to  the  start  of  the  observation  Interval  to  Improve 
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Che  esCloiate  of  In  practice,  this  is  not  an  option  open  to  us.  The 

minimal  problem  which  has  physical  meaning  is  therefore  that  of  jointly 
estimating  and  (>  (even  though  the  latter  is  ultimately  without  interest). 

Computation  of  the  remaining  elements  of  the  Fisher  information  matrix 
for  u)  and  proceeds  exactly  as  in  Eqs.  (7)-(13).  The  result  is 


J - 


1 M A/  . 

i 2 - t ] 

i-1  i 


I 1 ^ 2 2 

I -i/,  -h  ) 

I i-1  i 


1 M Ai  2 2 

1-1 


The  determinant  of  J is  easily  computed 


4 

1 

i 


M A. 

T I 

i-1  *'1 


N. 


(14) 


4 M M A ^ 

Det  J - ^ Z Z 

12  i-1  k-1  1 k 


Hence  Che  minimum  mean  square  error  fur  frequency  estimation  Is 

12 


(15) 


M A ' 

T r — ^ 
^ N 
1-1  1 


(16) 


This  quantity  is  independentof  t^,  as  one  would  expect. 

If  the  signal  to  noise  ratio  at  each  sensor  is  the  same 

. 2 


Then 


D^£o> 


2 

1 

12 


3 2 

T-*  M A /N- 


(17) 


M A^/N^  measures  the  output  signal  to  noise  ratio  of  a conventional  beam- 
former.  A single  sensor  with  that  signal  to  noise  ratio  would  there *ore 
be  entirely  equivalent  to  the  array  for  frequency  measuring  purposes. 


' • ' "v-  * .-rtA.*/.  1 . ’ 


r ^ 
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Next  we  must  look  into  the  coupling  between  frequency  estimation  on 
the  one  hand  and  bearing  and  range  estimation  on  the  other.  Since  bearing 
and  range  estimates  are  obtained  from  the  sensor  to  sensor  delay  of  the 
signal,  we  can  gain  the  required  Insight  by  working  with  only  two  sensors 
and  estimating  the  vector  At,  4>) , where  At  Is  the  signal  delay  between 

the  two  sensors. 

The  required  Fisher  Information  matrix  Is  readily  computed 


J - 


^ 


i“o  - h'j, 


2 ^ 
“o 


^ 

2 2^ 

- 'l  l I 

A I 


. if 

o 


^[(t  +T)2  - t 

Ni 


0)  T f 

o 


21^ 

“l 


(18) 


The  matrix  Inversion  Is  straightforward  and  one  obtains  the  mean 
square  estimation  errors 


2 * 6 

° <“o^  " -Ti 

T-*  A^/N, 


(19) 


2 A 2 

D^At)  - 


N, 


u ^ T 
o 


(20)  J 


Note  that  the  frequency  estimation  error  Improves  with  the  cube  of  the 
observation  ttee,  whereas  the  delay  estimation  error  only  Improves  with  the 
first  power  of  the  observation  time.  Perhaps  the  most  informative  compari- 
son is  that  between  the  fractional  errors  in  the  two  estimates 


This  ignores  the  ambiguity  problem  associated  with  sinusoidal  signals. 
See  Ref.  (!]. 
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9 “ 2 

D (u)  )/u 
o o 

D^(At)/(At 


(21) 


In  practice  T Is  likely  to  be  many  times  larger  than  At.  The  frequency 
estimate  therefore  becomes  exceedingly  accurate  before  a reasonable  delay 
estimate  can  be  made. 

Another  Interesting  comparison  Is  that  between  Eqs.  (19)  and  (17). 

Since  M - 2 In  Eq.  (19)  the  two  are.  In  fact.  Identical.  The  frequency 
estimate  Is  not  degraded  at  all  by  lack  of  a priori  knowledge  of  the  relative 
delay  (l.e.,  knowledge  of  source  location).  One  would  expect  this  relation- 
ship to  be  reciprocal  and  Indeed  one  finds  for  known  [using  the  last 
two  rows  and  columns  of  Eq.  (18) ] 


D^(At)|  --2-^ (22) 

u known  A 2 „ 
o t;—  0)  1 

“i  " 

This  Is  Identical  with  Eq.  (20),  confirming  the  Independence  of  the 
frequency  and  delay  estimates. 

2.  Narrowband  Gaussian  Signal.  TW  >>  1. 

In  many  cases  the  signal  observed  at  the  receiving  array  Is  not  a pure 
sinusoid.  There  Is  experimental  evidence  that  It  can  often  be  modelled  quite 
effectively  as  a narrowband  Gaussian  random  process.  The  unknown  parameter 
Is  now  Che  center  frequency  of  the  power  spectrum,  associated  with  this 
random  process. 

For  the  moment  we  consider  observation  times  T large  compared  with  the 
inverse  bandwidth  of  the  process  (TW  >>  1).  If  the  bandwidth  of  the  noise 
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Is  at  lease  equally  large,  the  most  convenient  fonn  of  the  data  vector 
uses  Fourier  coefficients.  For  TW  >>  1 the  Fourier  coefficients  associated 
with  different  frequencies  are  uncorrelated  [3].  Thus  the  covariance  matrix 
K of  the  data  vector  contains  many  zero  entries,  a feature  which  can  be 
exploited  by  arranging  the  data  vector  as  follows 


X ■ [c 


11  12’ 


C2i,C22.' 


•’'^2M ‘^nl’V’*‘'W 


(23) 


Cj^  Is  the  Fourier  coefficient  of  frequency  , measured  at  the  1^^  sensor. 
The  covariance  matrix  K now  assumes  the  block  diagonal  form 


(24) 


Kj  Is  the  spatial  covariance  matrix  of  Fourier  coefficients  at  frequency 

ui . . Its  elements  are  of  the  form  E{c..c*  }. 
j jl  jk 

The  basic  theory  of  parameter  estimation  from  a Gaussian  data  vector 
was  developed  In  Ref.  [4].  The  general  result,  for  unknown  parameter  a.  Is 


D^(i)  - {Tr[(K‘^ 


(25) 


Tr[  ] stands  for  the  trace  of  the  bracketed  matrix.  It  is  evident  from 

dK  2 

Eq.  (24)  that  (K  ^)  is  block  diagonal  and  Its  trace  Is  therefore 
simply  the  sum  of  the  traces  of  the  separate  blocks. 


(26) 


V .» 
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Kj  can  alao  be  simplified  a great  deal.  The  signal  and  noise  components 
of  X are  statistically  Independent.  For  an  array  of  moderate  dimensions  the 
same  signal  and  noise  spectra  are  received  at  each  sensor.  In  such  a 
case  one  can  write 

Kj  - S(Uj)P(a)j)  + N(Uj)Q(Wj)  (27) 

S((»)  and  N(u)  are  respectively  the  signal  and  noise  spectra  received  at  any 
one  sensor.  P(u)  Is  the  spatial  covariance  matrix  of  the  signal  at  frequency 
u),  normalized  so  that  the  diagonal  elements  are  unity.  Q(u)  Is  the  similarly 
normalized  spatial  covariance  matrix  of  the  noise.  The  parameter  the 
center  frequency  of  the  signal  spectrum.  Is  contained  only  In  the  scalar 
function  S(  ).  Eq.  (26)  therefore  assumes  the  rather  simple  form 


2 -1  2 -1 
D^(ui^)  - { Tr([S(aij)P(u)j)  + N(wj)Q(Uj)]  ^ P((i»j))  > 


(28) 


If  all  of  the  spectral  functions  are  essentially  constant  over  frequency 
Intervals  of  the  order  Aw  ~ 2'ir/T  [this  Is  essentially  the  large  TW  assumption] 
the  J sum  can  be  converted  Into  an  Integral 


-2,''  . 2it  , 
D (w^)  - — { 


Tr([S(w)P(w)  + N(w)Q(w)]"^  P(w))^dw} 

CtU) 

. o 


(29) 


To  proceed  further  we  assume  that  the  signal  wavefront  Is  coherent  over  the 
entire  array.  The  signal  received  at  one  sensor  Is  simply  a time  shifted 
version  of  the  one  received  at  any  ocher  sensor.  Under  these  circumstances 
the  P matrix  has  rank  1 and  can  be  written  as  an  outer  product  of  vectors 

P - V(w)  V*(u)  (30) 


••  ' 
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V(u))  is  the  "steering  vector",  the  vector  of  relative  delay  of  the  signal 

* 

at  different  sensors.  V (oj)  is  its  transpose-conjugate.  From  a standard 
matrix  identity 


IS(u))V(a))V*(a))  + N(a))Q(u,)] 


1 

N(m) 


S/m)  Q~^(M)Y(a))V*(oj)Q~^(a)) 


The  scalar 

G(u))  5 V*(u))Q"^(u))V(a)) 


(31) 


(32) 


is  recognized  as  the  array  gain  of  the  beamformer  optimally  matched  to  the 
given  noise  field. 

Substituting  Eq.  (31)  into  Eq.  (29)  and  temporarily  dropping  the 
aigximent  ui  for  greater  simplicity  of  notation  one  obtains 


n2/*  \ 2Tr 

D (w  ) ■ :;: 
o T 


o N o S 

1 + - G 


^TT 

T 


sc 

I 


(1  - )]^Trt(Q"^P)2]dco}“^ 

l.f  c 


(33) 


Furthermore 


Tr[(Q~^P)^]  - Tr[Q“^V  V^*q"^V  V*] 

- G TrlQ"^V  V*]  - G Tr(V*Q“2v)  - G^ 


(34) 


It  follows  that 


dS(to)) 
dw  2 

— ) do.}* 

S(u)) 


(35) 


Note  that  the  spectral  properties  of  signal  and  nolle  ere  quite  general, 
as  la  the  spatial  structure  of  the  noise. 
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If  Eq.  (35)  says  several  Interesting  things  about  the  estimation  of 

(1)  The  estimation  error  varies  as  T The  comparable  estimation  error 

-3 

for  a sinusoidal  signal  varied  as  T . For  the  narrowband  signals  It  Is 
no  longer  true  that  frequency  can  be  estimated  much  more  rapidly  than 
differential  delay  (bearing) . 

(2)  Array  geometry  and  spatial  structure  of  the  noise  field  enter  Eq.  (35) 

only  via  the  scalar  G(u) , the  optimal  array  gain.  output 

signal  to  noise  ratio  of  the  optimal  beamformer.  One  can  therefore  beam- 
form  prior  to  frequency  estimation  without  loss  in  performance  as  long  as 
one  uses  the  beam-former  structure  familiar  from  the  optimal  detection 
problem.  Alternatively  one  can  think  In  terms  of  processing  the  output  of 
a single  "equivalent”  sensor  whose  signal  to  noise  ratio  Is  ^^^G(u). 

(3)  Aside  from  the  output  signal  to  noise  ratio,  Eq.  (35)  depends  only  on 


lA 


Then  Eq.  (35)  becomes 


D^(w  ) - 
o T 


M F"  exp[- 


(u-w  )‘ 


’ (— 


(w-O) 

o 


° 1+  M ^ expl- 

o 2a 


w-w  2.  ,-l 
Ov  dw) 

2 ^ 

0 


If  Mjj—  <<  1 (i.e.  the  post-beamfonning  signal  to  noise  ratio  nowhere  reaches 
o 

unity)  the  integral  is  easily  evaluated  by  ignoring  the  second  term  in 


the  denominator.  Thus  for  low  signal  to  noise  ratio 


S 2 

2T(M  ^) 
"o 


If  >>  1 (high  signal  to  noise  ratio  in  the  signal  band)  an  appropriate 


approximation  is 


M — exp[- 
o 


(u)-u)  )‘ 


S (w-u  ) 

^ expl-  ] 

o 2o 


1 for  u-u)  1 < Aw 
I o'  — 


S (w-w  ) 

M — exp[ — ] for  I w-w  I > Aw 

N „ 2 O 

o 2o 


Aw  is  the  value  of  w-w  at  which  the  second  denominator  term  reaches 

' O' 


unity,  i.e.. 


- ^ a / log(M^) 


The  integration  of  Eq.  (38)  can  now  be  performed.  It  turns  out  that  the 
range  |w-w^|  > Awcontrlbutes very  little  to  the  result  and  one  obtains 


D^(w  ) ^ — 

o'  „ A- 


l/i  T [log(t^)]- 


The  principal  difference  between  Eqs.  (39)  and  (A2)  is  the  signal-to-nolse 
dependance.  VThan  the  post-beamformlng  signal  to  noise  ratio  is  low  the 


log(  ) is  the  logarithm  to  base  e. 
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the  mean  square  error  varies  as  the  Inverse  square  of  that  signal  to  noise 
ratio.  For  large  post-beamformlng  signal  to  noise  ratios  the  dependence 
becomes  quite  weak.  We  note  that  2Au>,  the  quantity  which  determines  that 
behavior,  Is  simply  the  frequency  band  over  which  Che  spectral  level  of 
the  signal  exceeds  that  of  the  noise.  Its  value  Is  therefore  quite  critically 
dependent  on  the  shape  of  the  signal  spectrum  and  no  general  significance 
should  be  attached  to  the  logarithmic  form  appearing  In  Eq.  (42). 

Simultaneous  estimation  of  frequency  and  bearing. 


We  must  now  discuss  the  Interdependence  between  the  center  frequency 
estimate  and  Che  location  estimate.  Since  we  are  not  considering  differen- 
tial Doppler  shifts,  we  are  presumably  dealing  with  a single  array  whose 
dimensions  are  small  compared  with  the  distance  to  Che  source.  In  chat 
case  Che  array  furnishes  little  useful  Information  concerning  range  but 
the  differential  delays  yield  a very  direct  measure  of  bearing.  We 
therefore  pose  the  problem  of  Jointly  estimating  center  frequency  and 
bearing  6. 

The  basic  theory  of  Joint  estimation  from  a Gaussian  data  vector  Is 
worked  out  In  Ref.  [4].  For  the  two  parameter  case  one  finds 


D2(i^) 

9 kno%m 

1 - P 


(43) 


where 


Trl(K-\^)^  Tr((K-^  ||: 
o 


P 


(44) 
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p measures  the  coupling  between  the  two  estimates  and  (1  - p)  is  the 


degradation  in  the  estimate  due  to  lack  of  prior  knowledge  of  6, 


The  first  term  in  the  denominator  of  Eq.  (44)  is  simply  [D  (u^) 


0 known 


and  is  therefore  available  immediately  from  Eq.  (35).  Similarly,  the  second 


term  in  the  denominator  is  [D  (0)  ] . It  is  derived  in  Ref.  [4] 

*^0  known 


(p.  96)  and  is  simply  restated  here  for  reference. 


Tr[(K-^  If)  ] - 


JL  1 , 


where 


6(u)  5 2{G(u))- 


dV*(u) 


dV(o)) 


-J6 


dV(u) 


|V  (a))Q  "(m)  I } 


This  rather  cumbersome  expression  will  turn  out  to  be  of  only  minor  impor- 


tance in  the  end. 


We  now  turn  to  the  numerator  of  Eq.  (44).  Decomposing  K into  single 


frequency  blocksK^  and  proceeding  exactly  as  In  Eqs.  (2S)-(33)  we  obtain 


Immediately  (once  again  suppressing  the  argument  Uj) 


N do)  , 

— ^ f 


To  compute  ^ we  recall  that  the  source  bearing  d affects  the 


spatial  covariance  matrix  P of  the  signal,  but  not  Its  spectrum  S.  Hence 

-1  HP 


Using  the  matrix  Identity  (31)  one  obtains  after  a few  steps  of  algebra 
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k 

ff 

t 


8K.. 


S 


K-^  Zi  - S N_  Q-lp,Q-l  dP 

J 36  1 + ^ ^ ^ de 

ir 


(49) 


Now  combining  Eqs.  (47)  and  (49): 


ii 

du) 


-1  !!S  -1  ^ 

K.  k/  -r4 

j J 30 

o 


it|o 


“ IQ-'p  - 


~ (q'^p)^]q'^ 


1 +|€ 


de 


(50) 


However 


G 


(oTV  - q“^v  v*q’^Y  V*  - gq"^p 


(51) 


Hence 


-1  -1 

aT  3e 


dS 

„ doj 

Q-V^f 

(1^|g) 


(52) 


Eq.  (44)  requires  us  to  compute  the  trace  of  this  quantity  and  sum 
over  the  frequency  Index  J: 


3K,  , 3K 


Trk-^  K-^  ^ 

30.^  “"j  30  „2 


dS 

du 


(1  + |g)^ 


Tr(Q‘^PQ“^  ^) 


(53) 


Now 


Tr(Q"W^  ^)  ■ Tr[Q'-‘-VVV^(^  V * + V ^ ) ] 


de 


de 


- Tr(V*Q"^V  V*g~^  1^)  + TrlQ'^V  V*  Q-^V  ^* ) 


(54)- 


t-r 


It  follows  that 


Tr(Q'S"^  £)  - G Tr(vV^  S Q"^V)  - G Tr(^  V*  q‘^V]  - G ||  —2 


(55)' 


Using  the  Invariance  of  the  trace  under  rotation  of  terms  in  its  argument. 
^Uslng  the  fact  chat  Q is  independent  of  6. 
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Finally,  substituting  Eq.  (33)  into  Eq.  (33)  and  evaluating  the  .j  sum  as 
an  Integral,  we  obtain 


Tr(K“^  K'^  ^)  - 
3u)  36^ 

o 


G(a)) 


dS  (m) 
du 

o 

S(w) 


dG(to)) 

de 


do) 


(56) 


Combining  Eqs.  (33  ),  (43)  and  (56)  one  can  now  write  down  a formal  expression 
for  the  coupling  factor  p of  Eq.  (44).  However,  the  information  of  primary 
interest  is  furnished  by  Eq.  (36)  alone.  Specifically^  the  integrand  of  (36) 
contains  the  factor  dG(u)/d6.  If,  for  example,  the  noise  is  spatially 
incoherent,  so  that  G(a))  ■ M,  dG/du  *■  0 and  hence  p *•  0.  In  this  important 
case  there  is  no  coupling  and  hence,  according  to  Eq.  (43),  no  degradation 
of  the  frequency  estimate  due  to  lack  of  prior  information  concerning  the 
bearing.^  More  generally,  the  coupling  between  frequency  and  bearing  esti- 
mation la  weak  as  long  as  the  array  gain  is  only  weakly  dependent  on  bearing. 
In  practice  this  is  apt  to  be  the  case  for  noise  fields  which  are  more  or 

leas  isotropic  or,  to  a lesser  extent,  for  any  noise  field  which  does  not 

2 

cuctain  strongly  directional  components. 


3.  Narrowband  Gaussian  signal.  TW  <<  1. 


Thus  far  we  have  only  considered  observation  times  large  coiiq>ared  with 
the  inverse  bandwidth  of  our  narrowband  Gaussian  siganl.  From  a practical 
point  of  view  this  is  undoubtedly  the  most  interesting  case.  However,  it 
is  cettainly  true  that  during  the  initial  moments  of  any  observation  TW 


1 2 * 

Since  the  expression  for  D (6)  is  analogous  to  Eq.  (43)  the  converse  is 
also  true:  Lack  of  prior  information  concerning  u does  not  degrade  the 

bearing  estimate. 

2 

Note  that  array  geometry  also  enters  into  G(w).  For  a linear  array  and 
an  isotropic  noise  field,  for  instance,  dG/du  ■ 0 for  a broadside  target 
but  becomes  large  within  one  or  two  beamwldths  of  endflre.  A spherically 
symmatrlcal  array  in  a similar  noise  field  would  yield  dG/du  0 for  all 
baarings . 
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is  small.  In  certain  situations  (e.g.  large  signal  to  noise  ratio)  we 
might  perhaps  be  fortunate  enough  to  obtain  the  desired  frequency  Informa- 
tion without  ever  reaching  the  large  TW  mode  of  operation.  A brief  treatment 
of  the  problem  therefore  appears  In  order. 

Before  beginning  this  discussion  we  must  be  aware  of  some  fundamental 
limitations.  A sample  function  of  a narrow^'cnd  Gaussian  process  looks 
essentially  like  a sinusoid  of  slowly  varying  frequency  and  amplitude.  For 
times  short  compared  with  the  Inverse  bandwidth  its  frequency  and  amplitude 
do  not  change  significantly.  All  our  estimator  can  do  is  to  establish  the 
present  frequency.  It  has  no  information  about  the  range  of  frequencies 
which  will  ultimately  be  traversed  and  can  therefore  not  establish  a 
center  frequency.  To  make  a center  frequency  measurement  one  requires  a 
statistically  significant  sample  of  the  random  process,  hence  an  observation 
time  satisfying  TW  >>  1.  The  short  term  measurement  is  meaningful  only  if 
we  are  not  concerned  with  the  frequency  fluctuations,  i.e.  if  we  do  not 
seek  an  accuracy  greater  than  the  bandwidth  of  the  narrowband  process. 

This  might  be  an  acceptable  condition  if  the  bandwidth  of  the  source  is 
extremely  narrow. 

Assuming  that  we  do  wish  to  make  the  short  time  frequency  estimate,  we 
are  faced  with  the  problem  of  estimating  the  parameter  of  a signal  s(t) 
given  by 

a(t)  ■ A 8in(u^  t - ^)  (57) 

This  Is  precisely  the  problem  posed  in  section  (1)  with  one  exception: 

The  anq>lltude  A is  an  additional  unknown  parameter.  One  must  therefore 
deal  with  the  3x3  Fisher  information  matrix  of  the  triple(a)^,  A). 
Computations  proceed  exactly  as  in  section  (1)  and  one  finds  (for  the 
same  signal  level  A and  noise  density  at  each  sensor  and  for  Tu^  '»>  1) . 
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iLA-  [(t  +T)^  - t 
3 J 


^ A [(t  +T)^  - t ^1 
2 1 1 ^ 


M A ._,2  ^ 2, 

^Kt^+T)  - t^  ] 


a2 


The  amplitude  estimate  Is  evidently  uncoupled  from  the  other  two  and  one 


obtains 


2 " 12 

D^u,  ) ^ 


This  Is  Identical  with  Eq.  (17).  It  Is  Interesting  to  observe  that  the  short 

3 

time  frequency  measurement  does,  Indeed,  Improve  with  T until  the  inherent 

randomness  of  the  process  asserts  Itself  (when  TW  ->  1) . It  Is  also  Important 

-2 

to  keep  In  mind  that  Eq.  (59)  depends  on  A . In  the  short  term  observation 
of  the  narrowband  process^  A Is  a random  variable.  We  may  be  lucky  and 
encounter  a large  value  of  A In  the  observation  Interval,  enabling  us  to  make 
an  accurate  short  time  frequency  estimate.  On  the  other  hand,  we  may  be 
unlucky  and  encounter  a period  of  very  low  signal  level  so  that  little  useful 
Information  can  be  extracted  and  we  are  forced  Into  the  use  of  larger 
observation  Intervals.^ 


The  parameter  A has  a Rayleigh  distribution,  hence  values  near  zero  are 
not  exceedingly  unlikely. 
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III.  Differential  Doppler  Estimation 

When  data  are  gathered  at  several  widely  separated  subarrays,  there 
may  be  differential  Doppler  shifts  between  the  various  subarray  outputs. 

These  provide  Information  concerning  source  heading  and  velocity  and  are 
therefore  of  considerable  practical  Importance.  As  In  the  case  of  single 
frequency  estimation  we  distinguish  between  sinusoidal  signals  and  narrow- 
band  Gaussian  signals  with  TW  >>  1.  Since  narrowband  Gaussian  signals  with 
TW  <<  1 are  effectively  sinusoids,  we  Ignore  that  category. 

In  much  of  the  development  concerning  frequency  measurement  without 
differential  Doppler  shift  we  allowed  the  array  to  consist  of  an  arbitrary 
number  of  elements  In  arbitrary  geometrical  configuration.  In  each  In- 
stance we  found  that  the  array  contributed  to  frequency  measurement  only 
through  enhancement  of  the  effective  signal  to  noise  ratio.  The  output  of 
the  optimal  beamformer  contained  all  available  Information  concerning  the 
frequency  or  center  frequency.  Beamformlng  prior  to  frequency  measurement 
Is  therefore  not  only  a practical  convenience,  but  It  Is  actually  the  best 
possible  procedure.  We  have  pointed  out  before  that  the  beamformer  output 
may  be  treated  as  the  output  of  a single  "equivalent"  sensor  with  Improved 
signal  to  noise  ratio.  From  here  on  we  shall  consider  each  subarray  as  such 
an  equivalent  sensor  and  shall  use  the  terms  "subarray"  and  "sensor" 
interchangeably . 

1.  SinuBoldal  signal 

Differential  Doppler  estimation  from  sinusoidal  signals  Is  basically  ^ 

the  same  problem  as  single  frequency  estimation.  Suppose  that  we  have  two 

1 

The  "optimal"  beamformer  la  that  which  maximizes  Che  output  signal  to 
noise  ratio. 
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subarrays  with  beamformer  outputs  x^(t)  and  x^Ct)  given  by  the  equations 

x^(t)  - A sin  [u)^(t-Tj^)  - + nj^(t)  (60) 

X2(t)  - A sin[(u^+Aa))(t-T2)  - (p^]  + n2(t)  (61) 

The  unknown  parameters  are  now  the  frequency  the  differential 
Doppler  shift  Aw,  and  the  two  phase  angles  and  1^2’ 

A few  words  of  explanation  may  be  in  order  concerning  the  need  to 
introduce  a second  unknown  phase  angle.  Suppose  we  had  chosen  ~ (^2 
as  in  part  II.  Except  for  small  differential  delays  • T2  [which  are 

discarded  in  steps  such  as  the  transition  from  (12)  to  (13)]  we  now  know 
that  the  signal  components  of  x^(t)  and  X2(t)  have  zero  crossings  which 
are  aligned  at  t ■ 0.  By  observing  their  relative  displacement  at  t^, 
the  actual  starting  point  of  the  observation  interval,  we  could  draw  strong 
inferences  concerning  Aw.  Indeed,  formal  computation  of  the  estimation 
error  yields 


♦l“^2' 


6_ 

(tj^+T)^ 


(62) 


The  formally  correct  but  practically  spurious  t^  dependence  is  precisely 
the  same  as  that  which  motivated  the  introduction  of  the  unknown  phase  angle 
^ in  section  II-l. 

Once  this  pitfall  is  properly  avoided  by  the  use  of  separate  and  ^2 
the  computation  of  the  Fisher  information  matrix  for  (w^.  Aw,  0^,  (^2) 
precede  exactly  as  in  Section  II-l.  The  result  is 
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(63) 


where 


After  a tedious  but  straightforward  matrix  Inversion  one  finds 

2A 


(64) 


D^(Aco) 


(65) 


and 


2 ' 12 
D («„)  - -y- 


(66) 


Comparing  Eq.  (66)  with  Eq.  (17)  [M  - 2]  we  note  that  the  form  Is  the  same, 
but  that  3db  of  performance  have  been  lost  through  Introduction  of  the 
second  unknown  phase  angle.  The  error  In  the  Au  estimate  Is  twice  as 
large  because  there  Is  only  one  differential  frequency  whereas  two  opportun- 
ities exist  to  measure  ui^.  We  also  note  In  passing  that  Eq.  (62)  becomes 
Identical  with  Eq.  (65)  If  one  chooses  tj^  ■ T/2,  the  worst  value  according  to 
Eq.  (62).^  We  note  further  that  computation  of  differential  Doppler  shift  from 
separate.  Independent  frequency  measurements  at  the  two  subarrays  leads  to  pre- 
cisely the  same  mean  square  error  as  that  given  by  Eq.  (66).  Coherent  processing 
of  the  s'lbarray  output  provides  no  advantage. 


A similar  observation  could  be  made  about  all  computations  in  section  Il-l 
In  which  the  time  origin  is  an  apparent  factor. 
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Fig.  2 


We  are  now  dealing  with  a 
source  radiating  a narrowband 
stationary  Gaussian  process  with 
spectrum  S(u),  symmetrical  about 
a center  frequency  Two 

receiving  subarrays  [Fig.  2] 
subtend  a sufficient  angle  B at 
the  source  so  that  the  Doppler 
shifts  of  the  received  signal 
components  are  not  Identical. 

We  are  concerned  only  with  ob- 
servation times  satisfying 
TW  >>  1 because  only  in  that 
case  are  we  dealing  with  a 


problem  significantly  different 
from  that  discussed  In  the  previous  section. 

Our  data  vector,  as  In  section  II-2,  Is  Gaussian.  The  basic  expression 
for  estimation  error  therefore  remains  Eq.  (25).  At  this  point,  however. 


the  strict  parallel  with  the  zero  differential  Doppler  problem  ends.  In 
section  II-2  we  were  able  to  block-dlagonallze  the  data  covariance  matrix  K 
by  working  with  Fourier  coefficients.  We  shall  still  find  It  useful  to 
work  with  Fourier  coefficients,  but  we  no  longer  achieve  the  Ingaediate 
block-dlagonallzatlon  because  differential  Doppler  shift  can  cause  correlation 
between  Fourier  coefficients  at  different  frequencies.  Our  first  task  Is 
therefore  computation  of  the  general  data  covariance  matrix  K. 
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The  data  vector  conalsts  of  the  Fourier  coefficients 


r 2 . , 

X,(w  ) - X,  (t)e"'^“n''  dt  1 - 1,2 

1 n j 1 


The  elements  of  the  covariance  matrix  are 


E{X^(u.^)X*(u.j^)}  - dt  do  R^j(r,a)e"^K'' 


where 


Rj^j(t,o)  - E{x^(t)Xj  (a) } 


1 - 1,2 


Is  the  temporal  correlation  of  the  received  time  functions.  E{  } stands 
for  the  expectation  of  the  bracketed  quantity. 

When  1" j , Eq.  (69)  Is  the  autocorrelation  function  of  x^(t)  and  Its 
computation  from  the  signal  and  noise  spectra  Is  a simple  matter.  We 
therefore  concentrate  on  the  case  1-1,  J-2. 

The  signal  coBq>onents  of  x^(t)  and  X2(t)  are  stationary  Gaussian 
processes  and  can  therefore  be  represented  by  series  of  the  form  [5] 


S,(t)  ■ I C cos  [(b)  + b)_>t  - ♦ ] 

1 q D'  ’^q' 


S2(t)  ■ E c^  cos[(b)^  + U)p  + Ab))(t  - T^)  - ♦q] 


The  are  statistically  Independent  random  variables,  each  uniformly 
distributed  over  (0,  2n).  c^  specifies  the  amplitude  of  the  sinusoid 
in  terms  of  the  known  signal  spectrum  S(b)) 


S(h)^) 


T 

o 
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is  a very  large  but  otherwise  arbitrary  interval  for  which  (70) 
and  (71)  are  defined.  For  our  purposes  we  require  only  >>  T.  The 
frequency  components  are  given  in  terms  of  (Tu  by  the  equation 

w ■ qrf'w  “ (73) 

Uq  is  the  Doppler  shift  at  subarray  1,  Au  the  differential  Doppler  shift 
between  the  two  subarrays  and  the  differencial  delay  of  Che  signal. 

The  representations  (70)  and  (71)  imply  a subtle  but  physically  im- 
portant assu]i9>tion:  Each  component  of  Xj^(t)  [and  similarly  X2(t)]  is 
Doppler  shifted  by  the  same  amount.  Actually  the  higher  frequency  com- 
ponents are  shifted  more,  thus  resulting  in  a slight  distortion  of  Che 
spectrum  which  is  being  Ignored.  The  magnitude  of  this  effect  is  bounded 
above  by  the  Doppler  shift  of  the  maximum  modulation  frequency  (signal 
bandwidth).  In  practice,  our  differential  Doppler  measurement  will 
depend  on  Che  coherence  of  s^(c)  and  a down-shifted  version  of  X2(c). 

That  coherence  will  be  affected  seriously  when  the  observation  time  T 
exceeds  a period  of  Che  differential  Doppler  shift  of  Che  highest  modulation 
frequency.  Thus  if  o is  Che  signal  bandwidth  in  rad/sec,  Av  the  difference 
of  Che  source  radial  velocities  measured  at  Che  two  subarrays,  and  c the 
velocity  of  sound,  our  computation  is  valid  for 

T — 0 « 1 (74) 

c 

For  values  of  T above  this  bound  the  notion  of  differential  Doppler  shift  is 
no  lunger  meaningful.  One  would  have  to  estimate  Av  directly,  implying  the 
use  of  true  time  compression  rather  than  mere  frequency  shift. 

Before  proceeding  with  the  computation  we  should  clarify  one  additional 
point,  in  section  lii-1  we  were  forced  to  introduce  the  random  phase 
*1*8'^*  $2  ovder  to  eliminate  the  spurious  dependence  of  Eq.  (62)  on  the 

,,■■■■  J I II  - 


a 
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time  origin.  It  Is  tempting,  but  Inappropriate^  to  Introduce  a similar 

random  phase  angle  Into  Eq.  (71).  To  do  so  would  mean  shifting  each 

frequency  component  of  (71)  by  the  same  angle.  Such  a procedure  changes 

the  waveshape  and  therefore  affects  the  correlation  between  (70)  and  (71). 

The  effect  becomes  serious  as  soon  as  TW  > 1,  precisely  the  region  In  which 

we  wish  to  carry  out  our  analysis.  This  leaves  us  with  the  problem  of  time 

origin  dependence.  We  shall  resolve  It  Initially  by  choosing  the  time 

T 

origin  at  the  center  of  the  observation  Interval  (t^  ■ -j) , the  choice 
which  yielded  the  correct  result  In  the  case  of  sinusoidal  signals.  Later 
we  shall  take  the  more  fundamental  approach  of  allowing  to  be  unknown, 
thus  eliminating  the  special  significance  of  the  time  origin.  We  shall 
find  that  Introduction  of  this  additional  unknown  parameter  does  not  change 
the  estimation  error  of  Am,  thus  confirming  the  validity  of  the  earlier 
result . 

In  much  of  the  subsequent  development  computational  details  become 
quite  cumbersome.  They  are  therefore  relegated  to  a series  of  appendices 
of  which  the  first.  Appendix  A,  deals  with  the  derivation  of  R^j(t,o)  from 
Eqs.  (69),  (70)  and  (71).  According  to  Eq.  (A-3) 


Rj^2(t,'7)  ■ |**  da)S(u))  cos [ (w+Wjj)  (t-o)  - Auo  + (w4-u)jj+Au)t^]  (75)  > 

0 

We  note  that  our  spectrum  S(ut)  is  so  defined  that 
«» 

1 S(w)du  ■ Average  signal  power  (76) 

o 

Only  the  signal  contributes  to  R^2  because  the  noises  at  the  two  subarrays  ^ 

are  assuaed  to  be  uncorrelated. 

; I 

Perhaps  the  most  striking  feature  of  Eq.  (75)  Is  the  fact  that  it  depends 
on  both  o and  t,  not  only  on  their  difference.  The  source  signal,  and  the  || 

signals  received  at  each  subarrey  ere  stationary  random  processes,  but  they  1 
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► * 


are  not  statlonarlly  correlated  because  of  the  differential  Doppler  shift. 
This  Is  at  the  core  of  the  difference  between  center  frequency  estimation 
and  differential  Doppler  estimation.  We  shall  find  that  the  periodic 
nature  of  the  non-statlonarlty  can  be  exploited  to  obtain  an  essentially 
coherent  estimate  of  Au. 

The  next  step  Is  the  computation  of  the  data  covariance  matrix  K from 
Eq.  (68).  Details  are  contained  In  Appendix  B.  Assuming  only  TW  » 1 
the  general  result  is  [from  Eqs.  (B~19)  (B-20)  and  (B-21)]: 


E{Xi(a)^)X2(w^)}  - 


irT 


sin(a).-u)  -Aa))';r  w.+u) 

Jl  n 2 I n 


Ab) 


s(- 


(‘*'r“n-A“)  2 


2 -“d>  ® 


, 'irT[S(a)  -w  ) + N,(u  )]  n - £ 
E{X,(w  )X*(a.,)}  - ^ ° ^ " 


E{X2(Wn)X2(V>  " 


f 7rT[S(a)^-Ujj-Au))  + ^2(1*)^)  n • £ 


1.  0 


n T*  £ 


The  matrix  K specified  by  Eqs.  (77)-(79)  Is  cumbersome  because  the 
cross-correlation  for  Fourier  coefficients  associated  with  different  sub- 
arrays [Eq.  (77)]  Is,  In  general,  different  from  zero  for  all  combinations 
of  n and  £.  Considerable  simplification  results  If  values  of  Aco  are 
confined  to  the  discrete  set 


Au  ■ , k an  Integer 


(77) 


(78) 


(79) 


Since,  by  definition 
_ 2TTn 


T ’ “£ 


2Tr£ 

T 


n,  £ Integers 


(80) 


(81) 


one  obtains 
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sin(aj.-w  -Aa))-;r 
X,  n £ 


Hence  Eq.  (77)  becomes 


8ln(t-  n-k)TT 

(a-n-k)TT 


0 


nT  e 


^ rnT  S (uj 


It  (u)  +Au)) 
-'o  n 


i - n+k 
I ^ n+k 


■ n+k 


1*  n-»k 


(82) 


(83) 


The  assumption  that  Au  takes  on  only  a discrete  set  of  values  2iTk/T 
should  be  of  small  practical  Importance.  Since  TV  >>  1,  the  allowed 
differential  Doppler  shifts  are  separated  by  much  less  than  the  signal 
bandwidth.  It  seems  unlikely  that  estimator  performance  for  this  rather 
dense  set  of  Au  values  should  not  be  generally  representative. 

For  later  manipulations  we  shall  find  It  convenient  to  arrange  the 
data  vector  as  follows 


(84)^ 


Then  the  form  of  the  covariance  matrix  K Is 


(85) 


The  complex  Fourier  representation  requires  both  positive  and  negative 
frequencies.  However,  negative  frequency 

coefficients  do  not  contribute  new  Information  and  can  therefore  be 
omitted. 


^ ^ ,v. 


**^^1  wi  ■v.jy  < k 
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j 

The  solid  diagonal  lines  Indicate  the  only  non-zero  entries  In  the  matrix. 

Coupling  between  the  two  subarray  outputs  exists  only  at  frequencies  separated 
by  the  differential  Doppler  shift. 

Because  of  the  sparse  form  of  Eq.  (85)  the  calculation  of  K ^ Is  not  | 

i 

difficult.  Computational  details  are  given  In  Appendix  C and  the  result 
Is  [for  spectrally  flat  noise  of  power  level  at  each  subarray] 
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where 


B - 

n 


D - 

n 


and 


S(u)  -<*»„)  + N 
n D o 

2 N Siui  -O  + N ^ 
o n D o 


2 N„  S 


p “ («  + Aw)  T 

n n o 


B.  Estimation  of  Differential  Doppler  Shift 

The  next  step  in  the  evaluation  of  £q.  (25)  requires  calculation 
of  the  derivative  (dK/da).  For  the  moment  we  shall  assume  that  Au  is  the 
only  unknown  parameter.  We  must  therefore  calculate  the  elements  of 
IdK/d(Aa))]. 

While  we  are  ultimately  concerned  only  with  the  discrete  set  of 

differential  Doppler  shifts  Aw  - 2iTk/T,  we  cannot  start  with  Eq.  (83) 

because  it  does  not  describe  the  Aw  dependence  near  the  point  of  Interest, 

which  is  required  for  the  computation  of  the  derivative.  We  must  therefore 

return  to  Eqs.  (77)-(79),  differentiate  with  respect  to  Aw  and  then 

evaluate  at  Aw  - 2iTk/T.  The  operation  is  trivial  for  Eqs.  (78)  and  (79). 

The  computation  for  Eq.  (77)  is  shown  in  Appendix  D.^ 

.2 


(87) 


(88) 


(89) 


1 j2Tr;ir- 

-Hj-?  S(w  -wj  - ^'(w  -wj]e 


d(Z^^*l<V*2<V^  - ' 


A-U 


T -'”n  “D 
i-n-k 


n “D' 


i-n-k 


_ (n+i-Hc) 

S[f(i+n-lO -wjj]e 


I ■ n+k 


i ^ n-tk 


(90) 


(91) 


32 


1-  nT  S ' (u)  -a)_-Au)) 
\ n D 


n - It 


n i<  H 


(92) 


The  various  components  of  K have  very  different  orders  of  magnitude. 

Since  « T the  dependent  term  In  the  first  line  of  Eq.  (90)  Is 
several  orders  or  magnitude  smaller  than  the  term  In  the  second  line.  The 
physical  Interpretation  of  (2Tr/T)S' (lo^-u^)  Is  the  change  In  S over  an  Interval 
of  (2Tr/T)  rad/sec  near  For  TW  >>  1 and  smoothly  varying  spectra 

this  Is  a very  small  fraction  of  the  value  of  S In  the  same  neighborhood. 

Both  terms  In  the  first  line  of  Eq.  (90)  are  therefore  orders  of  magnitude 
smaller  than  the  term  In  the  second  line  and  can  be  Ignored  without  Incurring 
significant  error.  The  only  non-zero  term  In  Eq.  (92)  Is  of  the  same  order 
as  the  S'  term  In  Eq.  (90).  It  will  therefore  be  Ignored  also.^  This 
leaves  the  following  simplified  version  of  (dK/d(Ab))). 


d(Au)) 


E{Xj^(u)^)X*(w^)}- 


1 - n-Ht 


(-1)^ 

2 l-n-k 


n-k 


(n+l+k) 

S[Y(l+n-k)-oajj]e  1 f n+k 


(93) 


(94) 


The  form  of  the  complete  dK/d(A(D)  matrix  is  shown  In  Eq.  (95) 


Mote  that  Eq.  (92)  characterizes  our  ability  to  estimate  Au  from  observa- 
tions at  the  second  subarray  alone.  This  Is  what  we  are  discarding 
when  we  make  Che  proposed  approximation. 
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(95) 


As  in  Eq.  (83),  the  solid  diagonal  lines  indicate  the  only  non-zero 
entries  in  the  matrix.  We  note  that  the  factor  (it-n-*k'  ^ in  Eq.  (93) 
causes  the  magnitude  of  these  non-zero  entries  to  vary  inversely  with  their 
distance  from  the  diagonal  line  of  zeros. 


— 1 H V ~ 

We  must  now  coi^ine  Eqs.  (86)  and  (95)  to  obtain  Tr{(K  dA(a))'^ 


Computational  details  are  shown  in  Appendix  E.  From  Eq.  (E-21) 


r5o  r^i 


- NQ[S(a)-  ^ -u)jj)  + S(a)  + ^ - a)jj)]  + 

Q [2NqS(u,  - ^ - V + Nq2][2NoS(u.  + ^ - V + 


S (u-Ujj)  do) 


(96) 


Eq.  (96)  can  be  used  for  a cuoibersome  but  rather  precise  computation 

of  the  Aui  estimation  error.  A much  simpler  form  can  be  obtained  by  making 

2 

an  additional  approximation:  Because  of  the  factor  1/r  only  small 
values  of  r contribute  significantly  to  the  sum.  For  small  r,  S(u)  will 
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not  vary  greatly  over  an  interval  of  jrr/T  and  one  can  approximate 


S(u)  + ^ - ujjj)  = S(u-Wjj).^  In  that  case 


■ 73  < \ 


S 


0 2No®‘"-“d>  "■  "0 


du 


The  mean  square  estimation  error  is  simply  the  inverse  of  this  quantity 

2 

The  r sum  of  Eq.  (97)  has  the  numerical  value  u /6  [6].  Hence 


D^(Aw) 


3 ® 


T 

12v 


S (ui-uq) 


2S^S(^-u>^)  + Nq 


dw 


C.  Comparison  of  differential  Doppler  estimation  with  center 
frequency  estimation  and  frequency  estimation  for  sinusoidal 
signal. 


I 


Perhaps  the  most  striking  difference  between  differential  Doppler 

estimation  [Eq.  (98)]  and  center  frequency  estimation  [Eq.(35)]i8  the  T 

2 A _3  2 

dependence  of  the  mean  square  error.  D (Au)  varies  as  T whereas  D (u^) 
varies  as  T In  this  respect  the  differential  Doppler  estimate  behaves 
like  the  frequency  estimate  of  a sinusoidal  signal  [Eq.(19)].  This  is  not 


The  approximation  is  better  than  one  might  think  at  first  glance.  If  the 
spectrum  varies  linearly  over  the  significant  range  of  tirr/T,  the  numerator 

2 

of  Eq.  (96)  Is  precisely  [ 2NqS (w-Up)  + Nq  ] . The  denominator  differs  from 

[2NQS(u)-mjj)  + Nq^]  only  by  4Nq^[S(w  - ^ - a)jj)S(a)  + ^ - ui^)  - (a.-u)jj)  ] . 

For  large  S/N  the  fractional  error  Is  therefore  oioportlonal  to  the  square 
of  the  change  In  the  spectral  level  divided  by  S^(w-uiq).  For  small  S/N  the 
approximation  la  even  better. 


- "‘I*’-*-  7 II 


(97) 


(98) 
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unreasonable:  Even  Chough  Che  signal  Is  a Gaussian  random  process,  che 

waveshapes  received  ac  Che  cwo  subarrays  are  deCermlnlsclcally  relaCed 

and  Che  differencial  Doppler  esclmace  can  therefore  proceed  on  a coherent 

basis. ^ No  such  possibility  exists  for  the  center  frequency  estimate  which 

Is  essentially  a measurement  of  power  distribution  over  frequency. 

To  explore  further  the  apparent  similarity  between  differential 

Doppler  estimation  and  frequency  estimation  for  a sinusoid  we  consider 

separately  the  cases  of  high  and  low  signal  to  noise  ratio  in  the  slgneil 

band.  From  a practical  point  of  view  the  high  S/N  case  Is  probably  the  most 

Interesting  one.  In  any  event,  It  Is  most  directly  comparable  with  the 

sinusoidal  frequency  estimation  problem  since  the  In-band  signal  to  noise 

2 

ratio  Is  almost  certainly  high  In  the  latter. 

a.  High  signal  to  noise  ratio  In  the  signal  band. 

When  S(u)  >>  N^  In  the  signal  band,  Eq.  (98)  can  be  approximated 
to  a good  degree  of  accuracy  by  the  equation 


D (Aw)  = 


24iiN, 


S(w)dw 


24Tr 


(99) 


Is  the  total  signal  power.  We  note  that  the  error  depends  only 
on  the  total  signal  power,  not  on  detailed  spectral  properties  of  the 
signal. 


word  of  caution  Is  In  order  here.  As  pointed  out  earlier,  we  have 
assumed  that  each  frequency  component  received  at  a subarray  Is  Doppler 
shifted  by  the  same  amount.  In  reality  the  higher  frequency  components 
are  shifted  by  slightly  larger  amounts.  This  causes  waveshape  distortion 
and  if  this  distortion  Is  not  the  same  at  both  subarrays,  there  will  be  de- 
correlation.  Time  compression  of  one  subarray  output  should  compensate 
for  tnis  effect,  but  the  analysis  presented  here  does  not  deal  with  that 
problem. 

2 

In  the  sinusoidal  case  the  effective  bandwidth  is  2n/T. 
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Eq.  (99)  is  not  yet  directly  comparable  with  Eq.  (17).  The  noise 
spectral  density  in  Eq.  (17)  was  defined  as  the  power  per  Herz, 
whereas  Nq  in  Eq.  (99)  is  the  power  per  radian/sec.  Thus 


N 


1 


2nN 


0 


The  sinusoidal  signal  power  in  Eq.  (17)  is 


With  this  nomenclature  Eq.  (17)  becomes 


12tt 

3 P 
T tfls 


(100) 


(101) 


(102) 


Since  the  two  sensors  (subarrays)  employed  in  differential  Doppler 
estimation  afford  only  one  opportunity  for  differential  frequency  measure- 
ment, one  should  probably  choose  M ••  1 in  Eq.  (102)  in  order  to  obtain  a 
fair  comparison.  With  that  understanding 

- 2 (103) 

“ <"o> 


Thus  the  differential  Doppler  measurement  for  a narrowband  Gaussian 
signal  differs  in  mean  square  error  from  the  frequency  measurement  for  a 
sinusoid  by  a fixed  factor  of  2,  Independent  of  signal  and  noise  properties. 

Another  and  perhaps  more  meaningful  comparison  can  be  made  between 
differential  Doppler  measurements  using  sinusoids  and  narrowband  Gaussian 
signals  respectively.  Prom  Eqs.  (65)  and  (99)  [using  Eqs.  (100)  and  (101)] 


D^(Au)  . . 

narrowband 

D^(Aw)  ^ 

sinusoid 


(104) 


37 


The  narrowband  differential  Doppler  estimator  is  therefore  equal  In  per- 
formance to  an  estimator  supplied  with  sinusoidal  signals. 

b.  Low  signal  to  noise  ratio  In  the  signal  band. 

When  the  signal  to  noise  ratio  does  not  exceed  unity  even  In  the 
signal  band,  Eq.  (98)  can  be  approximated  by 


(Am) 


127rN, 


S ((i)-Wjj)dw 


(105) 


In  contrast  with  the  high  signal  to  noise  ratio  situation  the  estima- 
tion error  now  depends  to  some  extent  on  the  shape  of  the  signal  spectnim. 
As  an  example  we  use  the  spectral  shape  of  Eq.  (37),  here  repeated  for 

reference  ~ 

(co-u)  )^ 

S(w)  ■ SQ{exp —2 } 


(106) 


A simple  computation  now  yields 


2 '' 
D (Aw) 


12i^ 


'low  S/N 


a ^0 


(107) 


N, 


To  convert  the  high  S/N  result  to  a comparable  form  we  note  that 

,2 


P - 

a 


f (w-Wq) 

S-  exp{ 5 — }du  ■ o Sq 

2o 


(108) 


Than  Eq.  (99)  becomes 


D^(A(ii)l 
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The  only  significant  difference  between  Eqs.  (107)  and  (109)  is  the 
signal  to  noise  ratio  dependence.  The  differential  Doppler  estimator 

_3 

operates  coherently  in  either  regime  [T  dependence] , but  its  performance 
figure  varies  as  the  square  of  the  signal  to  noise  ratio  for  low 
as  its  first  power  for  large  S^/N^. 

Finally  we  wish  to  make  a quantitative  comparison  between  the  differ- 
ential Doppler  estimator  and  an  estimator  which  obtains  the  differential 
Doppler  shift  from  separate  center  frequency  estimates  at  the  two  subarrays. 
Since  the  center  frequency  estimate  depends  on  the  signal  spectrum  we  continue 
to  work  with  the  example  of  Eq.  (106). 

a.  High  signal  to  noise  ratio  in  the  signal  band. 

From  Eq.  (42)  we  obtain  for  the  mean  square  error  of  the  center 
frequency  estimate  at  each  subarray  (M  - 1) 


3Tr 

iji  T 


o 


(log 


3/2 


(110)^ 


When  estimates  at  the  two  subarrays  are  subtracted  to  obtain  differential 
Doppler  shift,  this  mean  square  error  is  doubled.  The  direct  differential 
Doppler  estimate  yields  the  mean  square  error  given  by  Eq.  (109).  Hence 


D^ (Ai) 


D^(Aw) 


subtraction  of  estimates 


direct  estimate 


0.222 


(log  Sq/Nq) 


3/2 


(To)‘ 


Sq/Nq  » 1 


(111)^ 


^log  ( ) 


means  tha  logarithm  to  base  e. 
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The  most  important  feature  of  Eq.  (Ill)  Is  the  factor  (To)  . o Is 
the  signal  bandwidth  (in  rad/sec).  We  therefore  conclude  that  the  direct 
estimate  of  Au  is  better  than  an  estimate  derived  from  separate  center 
frequency  measurements  by  a factor  proportional  to  the  square  of  the  TW 
product.  A second,  though  lesser,  advantage  of  the  direct  measurement  Is  Its 
more  rapid  Improvement  with  signal  to  noise  ratio.  Since  the  logarithmic 
dependence  in  Eq.  (109)  comes  from  the  particular  spectral  form  of  Eq.  (106), 
the  signal  to  noise  dependence  could  vary  substantially  and  might  well  be  less 
favorable  to  the  direct  measurement  than  suggested  by  Eq.  (111). 

b.  Low  signal  to  noise  ratio  In  the  signal  band. 


When  Sq/Nq  1 ^he  mean  square  error  of  the  center  frequency  estimate 
Is  given  by  Eq.  (39).  For  M ■ 1 


.^2  / \ _ i/ii  0 

° sTi 


(112) 


The  mean  square  error  of  the  direct  differential  Doppler  estimate  Is  given 
by  Eq.  (107).  Hence 


(Aw) 


D^(Aw) 


subtraction  of  Wf.  estimates  „ 

II<Tor  . Sq/Nq  « 1 

direct  measurement 


(113)- 


Once  again  the  advantage  of  the  direct  measurement  procedure  Is  proportional 
to  the  square  of  the  TW  product.  For  <<  1 the  Improvement  Is  Inde- 

pendent of  the  signal  to  noise  ratio. 

Slote  that  TW  > 1 Implies  To  > 2it,  so  that  the  right  side  of  Eq.  (113)  is 
substantially  in  excess  of  unity  even  for  very  moderate  values  of  TU. 
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D.  Slmultaneoutt  estimation  of  differential  Doppler  shift  and  other  parameters. 


Thus  far  we  have  only  considered  differencial  Doppler  estimation  as  an 
Isolated  problem,  all  other  parameters  being  assiomed  known.  In  practice, 
this  Is  almost  certainly  unrealistic.  Parameters  such  as  source  bearing  and 
range,  signal  center  frequency  and  bandwidth  are  probably  not  known  a priori 
and  this  lack  of  knowledge  might  degrade  Che  differential  Doppler  estimate. 

The  amount  of  any  degradation  depends  on  Che  coupling  between  the  differential 
Doppler  estimate  and  estimates  of  the  ocher  potentially  unknown  parameters. 

We  now  study  this  coupling  on  a pairwise  basis. 

The  basic  relations  are  Eqs.  (43)  and  (44),  rewritten  here  for  Che 
estimation  of  Au. 


D^(Ai) 


a known 


1 - p 


(114) 


where 


t)] 


rxr(K-l  ^ k"^— ^ 

Sa  9 (Aw: 

„ ,,..-1  3K  ,,„-18K,^, 


(115) 


a Is  Che  second  unknown  parameter  whose  coupling  with  Aw  Is  under  Investi- 
gation. For  our  purposes  It  will  be  bearing,  range,  center  frequency  or  signal 
bandwidth.  Fortunately  It  turns  out  that  we  need  not  make  a separate  argument 
for  each  of  these  quantities. 

The  key  to  a simpler  argument  Is  the  manner  In  which  Che  second  para- 
meter a enters  Che  data  covariance  matrix  K whose  elemeits  are  given  by 
Eqs.  (77)-(79).  Bearing  and  range  affect  only  Che  delay  Tq.  Signal  center 
frequency  and  bandwidth  appear  only  In  Che  spectral  function  S(  ).  Differen- 
tiation with  respect  to  these  parameters  Is  a simple  procedure  which  does  not 


T 

f 
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affect  the  multiplying  factor 

T 

8in(u)j^-u)^-Au))-2 

(“rv^“)2 


In  Eq>  (77).  When  the  derivative  is  evaluated  at  Au  • 2irk./T,  k an  Integer, 
this  factor  causes  Che  differentiated  Eq.  (77)  to  vanish  unless  I * n-Hc 
Thus  9K/9a  retains  Che  form  of  Eq.  (85)  for  each  of  the  parameters  a under 
discussion.  It  is  now  a simple  matter  to  demonstrate  [see  Appendix  F for 
* X die 

details]  chat  K — also  has  the  same  form. 


Once  again  Che  solid  diagonal  lines  indicate  the  only  non-zero  entries  in 
the  matrix. 

3K 

The  ocher  factor  in  the  numerator  of  Eq.  (115)  is  K which  was 

previously  computed  in  Eq.  (E-9).  In  Chat  equation  B is  a diagonal  matrix  and 
A has  zeros  along  its  principal  diagonal  [see  Eq.  (E-8)].  Hence 


(116) 


I 


! 


■r#* 
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3 (Aw) 


has  the  form 


I 


(117) 


Comparing  Eqs.  (116)  and  (117)  we  note  that (117)  has  zeros  In  all  slots 
where  (116)  shows  non-zero  entries.  From  this  and  the  symmetry  of  the  two 
matrices  It  follows  Immediately  that 

If 

Thus  p - 0 and  from  Eq.  (114) 


D^(Ai)  - D^(Aw) 

a known 


Differential  Doppler  estimation  Is  therefore  completely  decoupled  from  the 
estimation  of  the  other  parameters  under  discussion.  It  Is  also  apparent  that 
this  conclusion  is  not  confined  to  the  four  specific  secondary  parameters 
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enumerated  above.  The  same  argument  would  apply  to  any  parameter  affecting 
only  element  to  element  delay  or  spectral  shape  of  the  signal.  Differential 
Doppler  shift  is  apparently  very  much  in  a class  by  itself  and  its  lack  of 
coupling  to  most  other  parameters  makes  it  a particularly  attractive  feature 
to  exploit  in  gathering  information  about  source  motion. 


Narrowband  Gaussian  signal.  TW  >>  1.  Three  subarrays. 


We  now  consider  a receiving 
array  composed  of  3 subarrays  which 
subtend  sufficiently  large  angles 
0 and  Y at  the  source  so  that  there 
can  be  measurable  differential 
Doppler  shifts.  We  deal  only  with 
the  case  of  narrowband  Gaussian 
signals  and  large  TW  products. 
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4^  Is  Che  differential  Doppler  shift  between  subarrays  1 and  2 and  Is  the 
time  delay  for  Che  same  pair.  ^2  is  Che  differential  Doppler  shift  between 
subarrays  1 and  3 and  T2  is  Che  corresponding  time  delay.  All  other  symbols 
retain  Che  definitions  established  in  connection  with  Eqs.  (70)  and  (71). 

Eqs.  (120) -(122)  also  retain  Che  earlier  assumption  that  differential  Doppler 
shifts  on  the  modulation  frequency  are  negligible. 

We  continue  to  represent  Che  received  data  by  Fourier  coefficients.  The 
data  covariance  matrix  K now  has  dimension  3N  x 3N  and  its  elements  are,  in 


cooq>lece  analogy  with  Eqs.  (77)-(79): 

E(Xi(u^)X2(«,))  - .T  S(  2 - Ujj)e 

at  n X ^ 


(123) 


E{X2(w^)X2(w^)}  - itT 
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(124) 
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E{X2(w^)X3(u,^)}  - nT ^ ^ ^ S(  -«n-^)e 


J-^(a.n+«-t+A2-Ai) 


(w^-u.n-A2+Ai)- 


^ rirT[S(u.  -a)_)  + N(«  )] 

E(Xi(Wn)Xi(V}  -j  o 


A ('TrT[S(u)  -(Dj^-A.)  +N(w  )] 

E{X2(u.„)X2(a.^)}  ^ 
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n « i 
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(125) 

(126) 

(127) 


E{X3*(a.^)X3(u.^)} 


f T[S(w^-Wjj-A2)  + N(u^)]  n - i 

I 0 n j*  i 


(128) 


As  in  ths  two  subarray  case,  drastic  simplification  can  be  achieved  by 
confining  the  actual  values  of  differencial  Doppler  shift  to  integral  multiples 
of  2v/T.  Thus  we  will  only  allow  the  actual  values 


46 


In  complute  anuiogy  with  thu  argument  leading  Lu  Kqu.  (93)  and  (94)  we 
obtain  [again  omitting  terms  of  order  T In  comparison  with  those  of  order 
T^] 


E{Xj^(Wjj)X*(w^)  } 


0 - n+k 


/ T (-1)^”"”^  IT  ir(n+l-lc) 


(135) 


ro 


E{X2(a)n)X*(U|^)  } - ■{ 


I " n+h-k 


/ „2  . ..i-n-h+k 

(_  L_  Izij ^.Ti 

^ 2 )i-.n-h+k 


^ ( a+n-h+k  ) -Wjj-A  ] e 


j^(n+£-h+k) 


i »*  n+h-k 


(136) 


Since  and  X^  do  not  contain  the  parameter 

^ M*1  0 <“T) 

Comparison  of  Eqs.  (135)-(137)  with  Eqs.  (131)-(133)  reveals  that  the  former 
are  zero  for  all  combinations  of  n and  I for  which  the  latter  are  different 
from  zero.  The  matrix  3K/3A^  therefore  has  zero  entries  In  all  positions  where 
Eq.  (134)  shows  non-zero  entries.^ 

Consider  next  the  derivatives  with  respect  to  the  delay  parameters. 

3K/3t^  and  clearly  have  the  same  form  as  Eq.  (134).  An  argument 

“1  9K  2 

analogous  to  Appendix  F shows  that  K - — , 1 ■ 1>2,  also  assumes  this  form. 

®^1 

^Ths  diagonal  blocks  of  3K/3A.  are  zero  to  the  order  of  our  approximation 
(sea  Eq.  (94)]. 

2 

Tills  la  almost  obvious  by  Inspection.  Since  K can  be  converted  Into  a block 
dlagooal  form  of  3x3  matrices  (Appendix  G) , K~^  and  3K/3t^  are  similarly  block 

dlfigonal  and  their  product  therefore  has  the  same  characteristic. 
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(138) 


exactly  as  in  the  case  of  two  subarrays.  The  lack  of  coupling  between 
differential  Doppler  estimates  and  estimates  of  most  other  potentially  unknown 
parameters  (bearing,  range,  signal  bandwidth,  etc.)  evidently  carries  over 
unchanged  from  the  2 subarray  problem.  In  discussing  differential  Doppler 
estimation  we  can  therefore  assume,  without  loss  of  generality,  that  the 
relative  delays  and  I2  ®re  known.  If  they  are  known,  we  can  certainly 
introduce  delay  elements  at  the  outputs  of  two  of  the  subarrays  which  align 
the  three  signal  components.  Since  these  relative  delays  are  reversible, 
their  introduction  cannot  alter  the  performance  of  the  optimal  differential 
Doppler  estimator.  In  evaluating  differential  Doppler  estimation  we  can 
therefore,  without  loss  of  generality,  set  ••  T2  “ 0,  thus  greatly  simply- 
fing  Eqs.  (131),  (132),  (135),  and  (136).  Calculation  of  the  mean  square 
estimation  error  can  now  proceed  in  straightforward  manner  much  as  in  the  2 
subarray  problems.  Computational  details  are  shown  in  Appendix  H.  With  all 
ocher  parameters  (including  A^)  known  and  with  the  same  approximations  as  in 
the  case  of  two  subarrays 
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(139) 


Eq.  (139)  is  somewhat  misleading  because  the  A^  and  A^  estimates  are  not 
uncoupled  and  the  assumption  that  A^  is  known  whereas  A^  is  not  appears  quite 
artificial.  We  must  therefore  study  the  degradation  of  the  A^  estimate  due 
to  lack  of  prior  knowledge  concerning  A2* 
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According  to  Eqs.  (114)  and  (115;  we  require  the  coupling  coefficient 


TrKK-l  ^)2]TrKK-'  ^)'] 


(140) 


The  first  factor  in  the  denominator  is  given  by  Eq.  (G-20) . Since  the  problem 
of  estimating  ^2  with  known  is  precisely  the  same  as  that  of  estimating  A^ 
with  A2  known,  the  second  factor  in  the  denominator  of  Eq.  (140)  is  also  given 
by  (G-20),  a fact  readily  confirmed  by  direct  computation.  There  remains  only 
the  calculation  of  the  numerator  of  Eq.  (140).  The  only  new  feature  in  it  is 
the  factor  dK/3A2.  The  elements  of  this  matrix  can  be  Inferred  immediately 
from  Eqs.  (135)-(137),  with  “ '^2  ” 
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(143) 


The  calculation  can  now  proceed  as  in  Appendix  G.  Details  are  given  in 
Appendix  H.  According  to  Eq.  (H-9)  the  coupling  coefficient  p has  a numerical 
value  of  (1/4)  completely  independent  of  the  spectral  properties  of  the 
signal.  The  mean  square  estimation  error  for  A^  with  A2  unknown  (or  A2 
with  A^  unknown)  is  now  easily  written  down  from  Eqs.  (139)  and  (H-9): 
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(144) 


It  is  Interesting  to  compare  Eq.  (144)  with  Eq.  (98),  the  latter  describing 

the  mean  square  estimation  error  of  the  differential  Doppler  shift  between  the 

signals  received  at  one  pair  of  subarrays.  When  the  signal  to  noise  ratio 

2 

in  the  signal  band  is  large,  so  that  the  term  in  the  integrand  can  be 
ignored,  Eqs.  (144)  and  (98)  are  identical.  Evidently  it  makes  no  difference 
whether  one  uses  the  subarrays  pairwise  to  obtain  separate  estimates  of 
and  A2  or  whether  one  works  with  all  three  subarray  outputs  simultaneously. 

When  the  signal  to  noise  ratio  is  low  even  in  the  signal  band  the  two  esti- 
mation errors  have  the  ratio 


D^(All)) 


1.5 


unknown 


(145) 


A modest  gain  can  therefore  be  made  in  principle  by  processing  the  three 
subarray  outputs  simultaneously.  It  appears  questionable,  however,  whether 
this  relatively  small  gain  would  Justify  the  increased  complexity  of 


instrumentation . 


Appendix  A.  Derivation  of  R^j(t,o). 


From  Eqs.  (69),  (70)  and  (71) 
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+ E{ C03  [ (w^+Wjj)  t-  (w^+iUjj+Auj)  0+  (lo^+ojjj+Aaj)  ■ro”‘*’q'^’^r 


(A-1) 


Since  rhe  ^'s  are  all  staClstlcally  independent  and  uniformly  distributed 
over  (0,2ir),  the  only  non-zero  contribution  is  made  by  the  second  term  in  (A-1) 
which  is  non-zero  for  q-r.  For  that  combination  the  (Ji's  vanish  and  the 
averaging  operation  becomes  trivial.  Hence, 

N c 2 

Rl2(tiO)  “ 2 ~2“  cos[(u)  +Wjj)(t-a)-Auo+(oj  +Wp+Aw)Tq]  (A-2) 

q-1  ^ 


Using  Eq.  (72)  and  recognizing  that  the  time  interval  Tq  can  be  made 
arbitrarily  large  (so  that  Aoi  0)  one  can  convert  the  sum  into  an  Integral 
00 

Rl2(t,o)  ■ J dw  S(ai)coB  [ (uj+Wjj)  (t-o)-Awa+((jjr*-u)jj+Au))  Tq]  (A-3) 

0 

In  completely  analogous  fashion 
00 

Rll(t,a)  ■ I da)[S(w)+Nj^(a»)  ]co8[  (uH-Wq)  (t-o)  ] (A-4) 

0 

00 

*22^*^’°^  ■ I dui(S(u))+N2(u))]co8[(u+u)Q+Aa))(t-o)]  (A-5) 

0 

Note  chat  the  autocorrelation  functions  (A-4)  and  (A-5)  incorporate  Che  received 
noise  spectra.  R^2  depends  only  on  the  signal  because  of  the  postulated  inde- 


pendence of  Che  noise  at  Che  two  subarrays. 
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Appendix  B.  The  elements  of  K. 


Substituting  Eq.  (75)  into  Eq.  (68)  and  using  the  exponential  form  of 
the  cosine  function: 


E(Xi(«^)X2(»j)) 


j (w+Ujj-o)^)  t f 2 -j  ((irhj)jj+Aw-u)^)a 


+ TT  du)S(u))e 
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(B-1) 


The  t and  o integrations  are  easily  performed 
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— J dwS(b))e 
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j (tD+u)jj+Au)TQ  sin(ui-fu)p-u)^)— 

X X 

(ur+Wjj-o)^)^  (urhi)jj+Au-u)j^)j 


a,  T T 

^2  -j  (urho^+Aw)TQ  sin(a)+Wjj+u)^)Y  8ln(uj+u)jj+Aw+w^)'2 

+ j-  du)S  (w)  e ^ 

■*0  ((i)+u)jj+u)^)-2  (u+Wjj-t-Aurhi)j^)-2 


(B-2) 


Our  signal  spectrum  S(u)  is  concentrated  in  the  neighborhood  of  u ■ Uq> 
The  second  integral  in  (B-2)  Is  therefore  essentially  zero  unless  Is  near 
Since  we  are  only  working  with  positive  frequency  Fourier  coefficients 


this  cannot  occur  and  we  can  concentrate  our  attention  on  Che  first  integral. 
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Define 
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Then  Eq.  (B-2)  becomes 
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dwe'^  *^S(u)(J)^j^(u)) 


(B-3) 


(B-4) 


We  have  extended  Che  lower  limit  of  integration  to  -<»  because  the  spectral 
function  S(b))  was  defined  as  zero  for  u)  < 0. 

Next,  using  F(  ) to  designate  the  Fourier  transform  of  the  bracketed 
quantity,  we  obtain  from  Parseval's  theorem 
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where 
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Also 
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The  ayvbol  iC  denotes  convolution. 
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We  know  that 


8in(urh.  -0,  )-  2n  T, 
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where  rect(Y;  j)  is  the  rectangle  function  defined  by 
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It  follows  that 


f T J(Y-p)(wD+Au-a)  ) 
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The  two  rectangle  functions  are  sketched  below  for  y > 0. 
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After  integration  and  a few  steps  of  algebra 
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(B-13) 
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Repeating  the  computation  for  Y < 0 one  finds  a change  only  in  the  sign  of 

y in  the  argument  of  the  sinusoid.  Thus  for  all  Y 

sinl(a,^-u.^-Aa,)|.(l  - M)]  j y ] 

*n^(y)  - e (B-14) 


Now  return  to  Eq.  (B-5) 


2 j(w  +Aa))T_  / 


dYRi(Y-TQ)«^^(Y) 


(B-15) 


R^(y)  is  closely  related  to  the  signal  autocorrelation.  In  fact,  suppose 
that  the  signal  spectrum  is  symmetrical  about  the  center  frequency  u ■ cuq  so 
that 

S((ju)  ••  G(ai-(.'Q),  G an  even  lunctlon.  (B-16) 


Then 


R^(Y) 


OP 

G(a>-u)Q)e 

^OD 


-JV| 


G(x) 


~*Yj 

e dx  ■ 


-Jv 


R(y) 


(B-17) 


where  R(y)  is  the  autocorrelation  of  the  signal  after  a down-shift  in  frequency 
by  u)q  rad/sec.  It  follows  that  R^  as  well  as  R vanish  for  arguments  in 
excess  of  the  signal  correlation  time.  For  TW  >>  1 this  time  is  much  smaller 
than  T.  Furthermore,  the  relative  delay  Tq  must  be  much  smaller  than  T if  the 
two  subarrays  are  to  process  information  coherently.  Under  these  conditions 
Y <<  T throughout  the  effective  range  of  integration  of  (B-15).  Hence  the 
sinusoidal  term  in  (B-14)  becomes  effectively  y Independent  and  one  can  write 


T 
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E{Xi(«n)X2(V^  - 


8in(u)j^-u)^-Au))-  j(Ujj+Aw)T^ 

T ^ 

(ta)  -w  -Aw)- 
A*  n z 


dYR^(Y-TQ)e 


-jY[“D  ■*■  2 ^ 


(B-18) 


The  integral  is  simply  a Fourier  transform  of  and  the  result  is  therefore 
easily  stated  in  terms  of  S(w).  After  a few  steps  of  algebra 


E{X^(w^)X2(Wj^)} 


T ^0 

8in(<D  -u)  -Aoj)»-  <*).+u)  -Aw  j— (w  +w  +Aw) 

ttT  = — S( w )e 

(w^-w^-Aw)2  2 


(B-19) 


The  remaining  terms  of  the  data  covariance  matrix  are  almost  obvious  by 
inspection.  When  both  Fourier  coefficients  come  from  the  same  subarray  output 
there  are  no  differential  delays  or  differential  Doppler  shifts.  The  signal 
spectrum  is  shifted  by  w^^  at  the  first  subarray,  by  w^  -t-  Aw  at  the  second. 
Correlation  now  exists  only  when  the  two  indices  are  the  same.  When  they  are 
the  same,  however,  there  is  correlation  for  the  noise  as  well  as  the  signal. 


The  correct  form  of  the  covariance  elements  is  now  evident  from  Eq.  (B-19). 


TrTlS(w„-w-)+N,  (w^)]  n-a 
n u in 


n j*  a 


(B-20) 


itTlS(w  -w„-Aw)+N«(w  )]  n ■ a 
n u z n 


n 1*  a 


(B-21) 


-A.. 
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Appendix  C.  Computation  of  K 


Rearrange  the  data  vector  as  follows 

X - »^2^“l+k^ ’ ’ " *^l^“N-k^ ’*2 I ^l^“N-k+l^ ' * " 1^*^^ I 

X2(u.^)...X2(u)j^)|  (C-1) 

In  this  arrangement  the  first  2(N-k)  terms  pair  all  indices  from  subarrays 
1 and  2 for  which  there  is  non-zero  correlation  according  to  Eq.  (83).  The 
remaining  terms  are  Fourier  coefficients  near  the  upper  edge  of  the  processed 
frequency  band  of  subarray  1 and  the  lower  edge  of  the  processed  frequency  band 
of  subarray  2.  Terms  in  the  latter  two  categories  are  uncorrelated  with  all 
other  elements  of  the  data  vector. 


The  data  covariance  matrix  now  assumed  the  block  diagonal  form 
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From  Eqs.  (78), (79)  and  (83)  the  2x2  matrices  are  given  by 


A ■ ttT 
n 


s(u)  -ui  )-m.  (u  ) 
n u In 


-J  Tq  (w^+A(i)) 


S(u)n-a>D)e 


jTo(w^+Aw) 


S(w  -6j_)+N_(u  +Au) 
n u z n 


(C-3) 


To  simplify  the  algebra  we  shall  assume  that  N^((ij)  N2((i))  ~ N^,  a 

constant.  Thus  we  are  postulating  that  the  noise  at  each  subarray  Is  spectrally 
flat  over  the  processed  frequency  band  wd  of  the  same  power  level  at 

each  subarray.  Since  the  signal  has  a very  narrow  bandwidth  and  any  Doppler 
shifts  will  be  small  compared  with  the  center  frequency  of  the  signal  this 

assumption  does  not  appear  to  Impose  any  serious  limitations. 


The  Inversion  of  A Is  now  trivial 
n 


A . 1 i 

ITT  2NqS(w^-u.jj)+N5 


S(<.^-a)p)-»-No 


-J  (u)^+Au))Tq 


-S(u>^-a)jj)e 


J (u)^+Au)Tq 


Stm^-a)jj)+No 


(C-4) 


The  matrices  B and  C are  diagonal.  The  frequencies  contained  In  B are 
all  within  Aw  of  the  upper  end  of  the  processed  band,  those  In  C within  Aw  of 
the  lower  end.  In  practice  one  would  certainly  process  a band  wide  enough  so 
that  no  significant  signal  components  are  lost  under  any  conceivable  Doppler 


■*T'' 
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shifts.  Hence  no  signal  power  should  lie  with  Au  of  the  band  edge  and  we  have 

B - C - ttTNqI  (C-5) 

where  1 Is  the  k x k Identity  matrix.  The  Inverse  of  Eq.  (C-2)  Is  therefore 


i(n-k)  2k 


[aV] 

as 

a;']  0 

0 

l( 

7 

0 - 
Fa’* 

- - r“-. 

1 

aTn,  0 

o 

• 

O * ' 

(C-6) 


Finally,  rearranging  the  data  vector  back  to  the  form  of  Eq.  (84)  [and 
hence  K to  the  form  of  Eq.  (85)]  we  obtain  the  result  shown  in  Eq.  (86). 
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Appendix  D. 


Computation  of 


dK 

d(Au))  ’ 


From  Eq.  (77) 


7m 


itT{- 


vT.T 


-(w.-o)  -Aa))7—  cos(u.-u)  -Au)-::+;r8in(u„-u)  -Arw)-:r  w.+w  -Aw  j-r— (w  +a).+Aa)) 

i— 2 2 L_Ei LJ i S(  ■.^■■■2  - -wje  ^ ^ 


(w.-w  -Aw)^ 

K n *7 — 


T 0 

sin(w  -u  -Au)-^  . u).+u  -Aw  t-  w +w  -Aw  j^w  +w  +Aw) 

+ I S-  w^)  + j-|  Sir^ Wj,) ]e  2 } 

(wrWn-Aw)2 


Now  write  Aw  ■ and  note  that  with  this  choice  (Wjj^-w^-Aw)-^  “ (i-n-k)ii^ 
It  follows  immediately  that  the  first  line  of  Eq.  (D-1)  vanishes  vdien  i-n-k  ■ 
whereas  the  second  line  vanishes  whenever  i-n-k  0.  Thus  Eq.  (D-1)  reduces 
to  the  relatively  simple  expression 


jT„(w  +Aw) 

f^tjTo  S(w^-Wd)-S' (w^-WD)]e  " i - n+k 


2 . .i-n-k  j-riv(n+i+k) 


(D-1) 


(D-2) 
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The  matrix  K ^ Is  given  by  Eq.  (86).  For  greater  ease  In  computation 
we  regard  It  as  composed  of  the  blocks  shown  In  Eq.  (E-1) . 


N-K 


^he  arrangement  Implies  that  k^  0.  This  Is  Insured  by  a convention  which 
applies  the  label  2 to  the  subarray  receiving  Che  signal  with  Che  higher 
center  frequency. 
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N-k 


JPl 


D e 
N-k 


■N-k 


(E-3) 


All  of  these  matrices  are  dlagonalf  B and  D of  dimension  (N—k^xCN—k) » N 
of  dimension  k x k.  The  blocks  labelled  0 in  Eq.  (E-1)  have  only  zero  entries. 

Next  we  represent  the  matrix  by  a similar  arrangement  of  blocks. 

Using  Eqs.  (93)-(95) 
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The  blocks  labelled  F,  G,  and  H have  Che  following  properties: 

1)  In  F and  H the  I Index  Is  between  1 and  k 

2)  In  G and  H Che  n Index  Is  between  N-k  and  N 


We  are  dealing  with  a narrowband  signal  concentrated  near  {See 

figure.]  The  Doppler  shift  and  differential  Doppler  shift  Au  > 2‘nk/T  are 
certainly  small  compared  with  Let  ~ 2:tR/T,  so  chat  R Is  the  Index  of 


\ 


of  the  Fourier  coefficient  associated  with  the  center  frequency  of  the  signal. 
Hence  the  factor  S[^(i+n-k)-Wjj]  in  Eq.  (93)  Is  different  from  zero  only  when 
Jl+n-k  is  near  2R. 

In  F and  H,  according  to  (1),  I Is  between  1 andk.  It  follows  that  n is 
of  the  order  2R  (compared  to  which  k Is  negligible) . Thus  Che  factor  i-n-k 
in  the  denominator  of  Eq.  (93)  has  a value  close  to  -2R.  From  the  definition 
of  R we  have 


2R 


(E-6) 


where  Tq  Is  the  period  of  the  center  frequency.  In  ell  interesting  situations 
the  ratio  T/T^  should  be  extremely  large  and  as  a consequence  all  elements 
of  F and  H are  asaantially  zero. 
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A very  similar  argument  applies  to  G.  Here,  according  to  (2),  n is  near 
N so  that  I must  be  near  2R-N  if  there  is  to  be  a non-zero  contribution. 

This  says  that  i-n-k  must  be  close  to  2R-2N.  N characterizes  the  upper 
edge  of  the  processed  frequency  band  and  will  certainly  be  very  much  larger 
than  R.  N can,  in  fact,  be  arbitrarily  large  so  that  the  elements  of  G 
are  arbitrarily  close  to  zero. 

Eq.  (E-5)  can  now  be  simplified  to  the  form 


dK 

d(Au)) 


0 ( 0 


0 t 0 


(E-7) 


The  elements  a of  A are  given  by  Eq.  (93)  with  i-k  = m.  Thus 
nm 


jTr^(n+ntf2k) 


m n 


(E-8)- 


1 2 

The  factor  T /2  in  Eq.  (93)  has  been  removed  from  the  definition  of  A by  its 
explicit  appearance  in  Eq.  (E-7). 
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Straightforward  matrix  multiplication  now  yields  from  Eqs.  (E-1)  and  (E-7) 


_1  dK 
^ d(Aa)) 


Therefore 


BA  \ N-k 


0 U k 


0 10  1 


-D*A  i N-k 


(E-9) 


^ d(Au))^ 


DA*DA* 

+ 

BABA* 


f-o-'t- 


BA*BA 

+ 

D*AD*A 


(E-10) 


Only  the  blocks  along  the  principal  diagonal  have  been  filled  in  since  we 
are  only  Interested  In  the  trace  of  the  matrix. 

The  required  trece  Is  now  easily  computed 


Xr((K“^  . JC 

d(Aw)^  ^ ^^2 


Tr(DA*DA*)  + Tr(BABA*)  + Tr(BA*BA)  + Tr(D*AD*A)} 


(E-11) 


Using  Eqs.  (E-3),  (E-8)  and  (89) 
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N-k  N-k  Jp  jp 

TrCDA^DA*)  - £ E D e D e a* 

n m mn 

n*l  m-1 


£ 

nq^n 


^0 

J2ir:jr-(n+nH-2k) 

D D e 
n m 


[^(iir+n)-a)jj 


-j2Tt“(n+nr»-2k) 

]e 


£ — - — 7 D D 

im^n  (m-n)^ 


S^[Y(m+n)-a)p] 


(E-12) 


Similarly 


Tr(D*AD*A)  - 


L £ 

n mf^n 


(m-n) 


2 


(E-13) 


Also  from  Eqs.  (E-2)  and  (E-8) 


Tr(BABA*)  - Tr(BA*BA) 


k-1  k~l 

£ £ B B a a*  - 

, , n m nm  nm 

n«l  in«l 


£ £ -VB„B  [^(n-Hn)  -o).  ] 

n m^tn  (m-n)^  n m T D 


I 


Hence 


Tr{ (k"^ 


t n m Ti  m T u 

mfn(m-n) 


(E-14) 


(E-15) 


(E-16) 
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Then  Eq.  (E-15)  assumes  the  form 


Tr{  (K 


-1 


d(AcL)) 


I — 
^ 2 


n r^O  r‘  [2NqS(u.^-<.jj)+Nq2]  [ZN^SCm^  + ^ + Nq‘] 


c2/  . T^r  V 

2,  5 ■*■  r -“d^ 


(E-18) 


The  n Index  formally  runs  from  1 to  N-k.  If  Che  arbitrary  upper  processing 
limit  N Is  taken  sufficiently  large  so  chat  all  signal  power  lies  below  ^(N-k) , 
the  n sum  can  be  extended  to  Infinity. 

For  large  TW  products  and  smoothly  varying  power  spectra,  changes 

only  Insignificantly  over  any  frequency  Increment  Am  » 2Tr/T.  The  n sum  can 
then  be  converted  Into  an  Integral 


3 ^ 2 

Att  r?*0 


Nq[S(u)*-u)jj)  + S(a>’  -1-2^  -oijj)]  + 


0 I2NqS(u.*-c0j,)  + No2][2NQS(a,*  + ^ -m^)  + 


S^(u)’  + ^ -(Ojj)  dm'  (E-19) 


2 

Because  of  the  rapidly  varying  factor  1/r  only  a few  terms  of  Che  r sum 
contribute  significantly  to  the  result. 

A somewhat  more  synmetrlcal  form  can  be  obtained  with  the  change  of 
variable 

u)'  + ^ - u)  (E-20) 


Since  there  Is  no  signal  power  near  w 0 Che  lower  limit  can  remain 


sero  and  Eq.  (E-19)  becomes 
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Appendix  F.  Computation  of  K 


Eq.  (E-1)  shows  K ^ decomposed  into  the  square  diagonal  matrices  B,D,N 

3K 

and  several  rectangular  null  matrices.  A similar  representation  of  — Is 

OOL 

according  to  Eq.  (85) 


^N-k 

U 

N-k 


3a 


The  matrices  P,Q,R,S,  and  T are  all  square  and  diagonal. 
Multiplication  of  Eqs.  (E-1)  and  (F-1)  yields 


v-1  M _ JL 

3a  itT 


-D*P+BR 


D*R+BQ 


^ N-k 


k 

k 

N-k 


(F-1) 


(F-2) 


Since  all  of  the  matrices  In  P-2  are  square  and  diagonal,  non-zero  elements 

occur  only  along  the  principal  diagonal  and  the  same  pair  of  secondary  diagonals 
—1  dK 

as  In  K..  Thus  K has  the  form  shown  In  Eq.  (116). 
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•"1  3K  2 

Appendix  G.  Computation  of  Tr{(K  } for  3 dubarrays. 


Rearrange  the  data  vector  as  follows  (let  h ^ k ^ 0) 

X - **3^“l+h^  ‘ ‘ '^l^“N-h^  ’*2^“N-h+k^  •h 

• i 

*l^‘^-k+l)  ’*2^‘^-h+k+l)  ’ • • ’^l^-k)  ’ • *l^“N-k+l^  ’ * * '*1%^  ’ | 


^2  (u)^)  » X^  (^2^  I ^ f • • • X2  (u)^)  f 


^3^***1^  I • • ' 


k) 


(G-1) 


All  frequencies  In  the  second  and  third  lines  of  this  expression  are  within 
A2  of  either  zero  or  As  In  Appendix  C we  shall  assume  that  the  signal  band 

Is  sufficiently  narrow  and  has  been  chosen  sufficiently  large  so  that  there 
Is  no  appreciable  signal  power  In  these  border  regions.  Then  the  last  two 
lines  contain  only  noise  (uncorrelated  from  sensor  to  sensor)  and  the  data 
covariance  assumes  the  rather  simple  form 
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The  M are  3x3  matrices  readily  obtainable  from  Eqs.  (126) -(128)  and 
n 


(131)-(133)  {with 


M - nT 
n 


S(u.^-a)jj)+No  I S(a)^-a)jj) 


S (w  -w  ) 


s(vV‘^o 


(G-3) 


The  Inverse  of  Is  easily  computed 


2S(V“d)-«'0 


”n^  ” ■itTNq(3S(u)^-u)jj)+NqJ 


-S(u)^-u,d) 


-S(u)^-u)jj)  I -S(VV 


2S(a.n-u)^)-WQ  j -S(w^-u>^) 

I 

4. 


-S(a)n-Wp)  I 2S(u)^-u»jj)+Nq 


(G-4) 


Now  rearranging  the  data  vector  back  to  Its  original  form  we  obtain  the 
following  explicit  version  for  the  Inverse  of  Eq.  (134). 
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Here 


(G-6) 


(G-7) 


All  entries  In  Eq.  (G-5)  other  than  the  Indicated  diagonals  are  zero. 

The  form  of  -rr-  Is  obtainable  directly  from  Eqs.  (135)-<137).  As  In 

9dl 

the  two  subarray  case  we  Ignore  the  "edge  effects"  near  the  upper  and  lower 
end  of  the  processed  frequency  band  [see  transition  from  Eq.  (E-5)  to  (E-7)} . 
Then 


K 


(G-8) 


C Is  an  (N-k)  X (N-k)  matrix  with  the  elements 


c 

nm 


m"n 


m ft  n 


D Is  an  (N-h)  X (N-h)  matrix  with  the  elements 

0 m«n 


nm 


S[^(m+n)-u)  ] m ft  n 
m-n  1 u 


Eq.  (G-5)  can  be  put  Into  a form  similar  to  (G-8)  with  the 


definitions 
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K K 


: 


-laK  _T 
3A^~2ir 


! 


C' 


1 

1 

' 1 -V 

A^“'  C 1 - 
1 ^ 

0 0 

1 

D"  1 

1 

1 

4 

-E?‘^D 


.__C ^___  o__; o_ 

O I o , o 


(H)  r ^ 


O 


I 

I 

0|0 

I 

I 

1-^-  -I 

I o 
I 


I 


ojo 


(k>^a) 

D C. 

+f^'^  D* 


O 


k 

If< 


D 


O 


O 


u 


-b‘% 


Jj  U-K 


'} 


U K 


k>K 


In  analogy  with  Eq.  (G-11) , Is  used  to  designate  the  matrix  C with 

the  last  h-4i  rows  and  columns  deleted . Tha  second  colximn  of  Eq . (G-13)  has 
componants  of  dlffarant  dlmanslonallty.  This  Is  to  be  Interpreted  as  follows: 
Before  addition  to  the  matrix  la  raised  to  dimension 

(M-k)  X (N*k)  by  adding  rows  and  columns  of  zeros  at  the  bottom  and  at  the 
right. 


(G-13; 


k 
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i 

i 


A straightforward  computation  now  yields 


, , -1  3K  ,2, 
Tr{K  — ) } - 
3Ai 


4it^ 


+ 2A^*'^DA^^^D*  + B^^^DB^^^D  + B^^^D*B^^^D*} 


(G-U) 


Apparently  Inconsistent  matrix  products  such  as  (B^^^D)A^^^C*  are  made 
meaningful  by  adding  rows  and  columns  of  zeros  at  the  bottom  and  right  of 
the  (N-h)  X (N-h)  matrix. 

The  problem  of  different  dimensionality  In  different  components  of  (G-14) 
becomes  trivial  when  the  trace  Is  written  out  as  a sum.  When  the  Indices  n 
and  n are  very  different  Eqs.  (G-9)  and  (G-10)  Indicate  that  the  terms  In 
question  are  small.  When  both  m and  n are  within  h of  the  upper  bound  N 
they  are  also  small  because  there  Is  no  signal  power  near  the  edge  of  the 
processed  band.  Thus  both  Indices  can  be  allowed  to  run  over  all  Integers 
and  Eq.  (G-14)  becomes 


— 1 iiK  7 

Tr{(K  ^ ~)  } - ^ ZE{2BBcc  +AAcc  -2BAdc 

3A,  ,2  nmmnnm  nmnmnm  nmmnnm 

1 4ir  n m 


>2BBcd  -«-2BBdc  -2ABdc  -t-ZAAdd 
nmnmnm  nmnmnm  nmnmmn  nmnmnm 


-«*  2B  B d d } 

n m nm  mn 


(G-15) 


According  to  Eqs.  (G-9)  and  (C-10) 


• -c 

mn  nm 


d ■ -d 


c 


(G-16) 
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hence 


Tr{K“^  — )^}  - ^ I U2A  A - 4B  B -BA  - A B }c^ 

3A'  -2  nm  nm  nm  nmnm 

1 2iT  n m 

Now  substituting  the  definitions  (G-6),(G-7)  and  (G-9)  and  making  the 
change  of  variables  m - n+r 


(0-17) 


iL  . i_  I s2[i(2n-h:)-«  1 


(G-18) 


This  is  the  equivalent  of  Eq.  (E-18)  in  the  two  subarray  problem.  To  simplify 
the  result  we  assume  again  that  the  spectrum  does  not  change  a great  deal  over 
the  significant  r range.  Then 


Tr{(K-^ 


•O  1 

( I E 
r“l  r n 


S (u)  — to  ) 
' n D' 


CG-19) 


FlT4ally,  converting  the  n sum  into  an  integral  and  inserting  the  numerical 
2 

value  TT  /6  for  the  r sum 


^ 3Kv2,  T 

Tr{(K  3^)  } - ^ 


S (w-tOjj) 


NQ[3S(a)-u)jj)+NQ] 


du 


(G-20) 


The  reciprocal  of  this  quantity  is  the  mean  square  error  given  by  Eq.  (139). 
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—1  —1  3K 

Appendix  H.  Computation  of  Tr(K  K 


3K 


We  begin  by  writing  — — In  the  form  of  Eq.  (G-8) . From  Eqs.  (141)-(143) 

dA/> 


^ 7j 
Z 

‘'{I 


r^i  -I 


o 

I + 

o 


! o 


c 


c 


--t^3 

0 o C 


4 o’ 

0 I 
I 


u 


Vx 

K 


^ K 


(H-1) 


The  elements  of  the  (N-h)  X (N-h)  matrix  C are  given  by 
'0 


nm 


m ■ n 


vin—n 


- — S[Y(m+n)-(A)-]  m »*  n 

iii*n  1 u 


(H-2)^ 


^In  terae  of  the  nomenclature  of  Appendix  G,  the  C of  (H-1)  is  of  G-13. 


From  Eqs.  (G-13)  and  (H-3) 
„2 


— 1 TkK.  —1  T* 
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From  Eqs.  (G-10)  and  (H-2)  it  is  apparent  that 


c “-c  ■ d ■ -d 

nm  mn  mn  nm 


(H-5) 


Hence  Eq.  (H-4)  becomes  (after  a few  steps  of  algebra) 


Tr(K“^  ^ E Z (4B„B  + A B + B A - 2A  A )c^ 

34,  dA«  ,2  n m n m n m n m nm 
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Now  substituting  the  definitions  (G-6),  (G-7)  and  (H-2)  and  making  the  change 
of  variables  m ■ n+r 
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(H-7) 


—1  3K  2 

This  expression  is  exactly  half  of  Tr{K  *9^)  ^ given  by  Eq.  (G-18).  Since 


Tr{(K  — ) ) . lr{(K  jj-)  ) 


(H-8) 


The  coupling  coefficient  p defined  by  Eq.  (140)  assumes  the  numerical  value 


(H-9) 
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